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Introduction to drive mechanisms 
Field development project 



► Field development projects begin with a phase consisting in acquiring, 
gathering and analyzing geological, petrophysical and fluid data 

► These data are necessary in order to: 

• Describe the reservoir as accurately as possible 

• Architecture: top, bottom, lateral extension and closure, any possible compartmentalization 
(faults), presence and distribution of the heterogeneities 

• Fluids: fluid properties, contacts 

• Estimate volumes in place/accumulations 

► Then drive mechanisms must be determined as early as possible in order to 
identify the way the reservoir is produced and how it will behave in the 
future 

• Estimate the recovery factor and the reserves 

• Determine the production scheme and corresponding production profile 

► Identifying the drive mechanisms is mandatory to help to make a decision for 
the Field Development Plan (FDP) 
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Introduction to drive mechanisms 
Definitions 



► Drive mechanisms 

• Mechanisms of production of the reservoir 

• How is the energy necessary to produce/recover the hydrocarbons provided? 


► Three types of drive mechanisms/recovery 

• Natural drive/Primary recovery: the field is produced through its own energy 

• Immiscible fluid injection/Secondary recovery: energy is provided to the field through 
injection: water flooding, gas flooding, cycling 

• Enhanced oil recovery/Tertiary recovery: energy for production is provided through 
complex techniques: 

• Miscible fluid injection 

• Chemical process: polymer flooding, surfactant flooding 

• Thermal process: steam flooding, in-situ combustion 
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Introduction to drive mechanisms 



• Monophasic expansion, compaction 
•Solution gas drive, gas cap drive 


•Water influx 
•Gravity drainage 



TERTIARY RECOVERY 


Enhanced Oil Recovery 



• Steam injection 

• In situ combustion 

• Electric heating 


• Miscible Hydrocarbons 

• C0 2 

• N 2 


• Polymer 

• Surfactant 

• ASP 
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Introduction to drive mechanisms 


\ 

► Drive mechanisms refer to the way the energy necessary to produce 

hydrocarbons is provided to the reservoir 

► Three types of drive mechanisms 

• Natural drive/Primary recovery: the reservoir is produced through its own energy 

• Immiscible fluid injection/Secondary recovery: energy is provided through fluid 
injection: water injection, gas injection, cycling 

• Enhanced Oil Recovery/Tertiary recovery: energy is provided through complex 
techniques 

► Identifying the drive mechanisms is mandatory in order to help make a 

decision for the Field Development Plan 

• Generally, several drive mechanisms may be taken into account, spanning over the 
whole life of the reservoir 

• IOR/EOR are generally considered from the beginning of the process to establish the 
FDP but identifying the natural drive mechanisms of the reservoir is always important in 
order to understand its behavior 

V I 
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Introduction to primary recovery 
Definition 



► Primary recovery refers to the recovery through the natural energy of the 
reservoir system i.e. natural drive mechanisms 

► There are various natural drive mechanisms coming from the impact of a 
number of physical phenomena: 

• Expansion of reservoir fluids: oil, gas and water 

• Liberation and expansion of solution gas (i.e. lighter part of hydrocarbons mixture) in the 
case of oil reservoirs 

• Expansion of reservoir rock and reduction of pore volume 

• Pressure support and sweeping from an adjacent active aquifer 

• Action of the gravity forces 


► Natural drive mechanisms have to be understood and evaluated as early as 
possible in the field history to assess the recovery factor and make a decision 
about the production scheme of the reservoir 
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Introduction to primary recovery 

Reminders: accumulation, reserves and recovery factor 



► Accumulation: volume of oil or gas initially in place 


► Reserves: recoverable volume of oil or gas 
Reserves = Accumulation xRf% 


► Reserves depend on a geological model, a reservoir model, a development 
scenario, economics, laws and contracts 


► All volumes of accumulations (OOIP, OGIP or GIIP) and reserves are expressed 

in surface (standard) conditions 

• Standard conditions: 1 atmosphere, 60°F 
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Introduction to primary recovery 


Reminders: phase envelope of a hydrocarbon mixture 
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Introduction to primary recovery 
Reminders: PVT parameters against pressure 



Evolution of the main reservoir properties with pressure 



Introduction to primary recovery 



Reminders: compressibility 

► Compressibility is related to the relative change in volume of an element (gas, 
liquid, solid) when it is submitted to a variation of pressure 


1 dV 
V dP 


dV = -cVdP 


► Typical compressibility values in the reservoir 

• Oil: c 0 = 0.7 to 3 x 10' 4 bar * 1 (above bubble point) 

• Water: c w = 0.4 to 0.7 x 10' 4 bar 1 

• Pore volume: c p = 0.3 to 1.5 x 10‘ 4 bar 1 

• Gas: c g = 7 to 145 x 10 4 bar 1 (highly dependent on the pressure of the reservoir) 

► Compressibility is the physical parameter behind fluid and pore volume 
expansion and compaction and one of the main parameters involved in 
natural drive mechanisms 
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Introduction to primary recovery 
Reminders: production parameters 



► Typical parameters of interest to monitor production and assess the 
performance of a recovery mechanism and related production scheme 


• GOR = Gas Oil Ratio 

• Instantaneous GOR: 


GOR = 


Gas flowrate@S.C 
Oil flowrate@S.C. 


• Average GOR: 


GOR av 



Produced Gas volume@S.C 
Produced Oil volume@S.C. 


Gp 

N p 


• There might be large differences between the two at a given moment of the production! 

• Average GOR is used for material balance calculations (see below) 


• WC = Water-Cut 

• Instantaneous WC: 


WC = 


Water flowrate@S.C 
Water flmowr ate + Oil flowrate@S.C. 
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Introduction to primary recovery 
Primary recovery for oil and gas reservoirs 

► Oil reservoirs show four main natural drive mechanisms that are split into 
two groups 

• Depletion drive: the reservoir is not in contact with a large active aquifer 

* Monophasic fluid expansion & pore volume compaction 

* Solution gas drive 

* Gas cap drive 

• Natural water drive: the reservoir is in contact with a large active aquifer providing 

pressure support and sweeping of the reservoir 

► Gas reservoirs are mainly concerned with monophasic fluid expansion 

• Depletion drive: very good performance 

• Natural water drive: poor performance! 

• Special case: gas condensate reservoir 

* Try to keep above dew point pressure! 
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General principle and performances 


► At pressures above bubble point pressure, oil, connate water and rock are the 
only materials present in the reservoir 

► As the pressure drops due to the production, the reservoir fluids expand 
while the pore volume compacts due to their individual compressibilities 

► As a result of expansion of fluids and reduction of pore volume, more oil (and 
possibly water) is forced out to the wellbore 

► Typical performances 

• Reservoir pressure declines rapidly while GOR remains constant and equal to R s , the 
solution GOR 

• Water-Cut is generally nil since this drive mechanism occurs very early in the life of the 
reservoir, when production is still anhydrous 

• Volume variations that are behind this mechanism are related to the compressibilities of 
oil, water and pore volume, which are generally quite low: as a consequence the 
expected recovery factor is very low, in the range of a few % only 
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Compaction drive 

General principle and performances 



► During natural depletion (no or weak pressure support), the fluid pressure FP 
decreases, hence the grain pressure GP increases and the pore volume 

shrinks according to pore compressibility: 

1 A(PV) 

Cp ~PV AP 


► Generally, c p is low (0.3 to 1.5 x 10' 4 bar _1 ) and assumed constant all along 
the monophasic natural depletion 

• Pore volume compressibility may be neglected for sake of simplicity leading to an 
assumption of constant pore volume 


► However, in some cases, typically reservoirs with high initial porosity, c p may 
increase leading to significant additional oil recovery 

• Compaction drive mechanism may result in Recovery Factor up to 20% 
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Compaction drive 

General principle and performances - 2 



► Pore volume compressibility and bulk rock compressibility 

• Bulk rock compressibility is the one that is measured in laboratory (through SCAL) and it 
is defined as: 

1 A (BV) 1 Ah 

Cb ~BV AP ~ h AP 

• As A (BV) = A(PV) we get: 

BV 

Cp Cjj ~ — Cjj/ 0 

► Generally high c p are associated to overpressure and/or high initial porosity 

• When fluids are produced out of the pore volume, the fluid pressure decreases leading 
to an increase in the Grain Pressure and a decrease in porosity thus forcing fluids out of 
the pore volume 

• The process generally results in a mechanical failure of the rock and a collapse of the 
porous matrix => compaction drive is often associated with surface subsidence 
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Compaction drive 
Pressures in the reservoir 



► Overburden Pressure 

• Overburden pressure is the sum of the fluid 
pressure and grain pressure 

OP = FP + GP 


Op Weight ( fluid + rock matrix ) 

Area 

• OP gradient is about 2.5 bar/m or 1 psi/ft 


Pressure (psia) 



► Fluid pressure 

• Fluid pressure or formation pore pressure is 
the pressure acting upon the fluids (water, 
oil, gas) trapped within the pores of the 
formation 

• FP gradient is about 0.1 bar/m or 0.45 psi/ft 
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Compaction drive 
Pressure regimes in the reservoir 

► There are three types of fluid pressure regimes 

• Hydrostatic pressure: the pressure corresponding to the weight of a vertical column of 
water up to surface 

• Subnormal pressure or underpressure: a fluid pressure less than the normal hydrostatic 
pressure 

• Overpressure: a fluid pressure superior to normal hydrostatic pressure 

► Overpressure 

• may be caused by compaction, formation uplifting, faulting, repressurization, 
hydrocarbon generation, massive evaporite deposition, mineral phase change etc. 
actually a rapid change in the strain system so that fluid pressure has not enough time to 
reequilibrate => generally associated to low permeability or permeability barriers 

• needs the reservoir to be sealed off from the surrounding strata => several pressure 
regimes may exist if reservoir is compartmentalized 

► Overpressure plays an important role in the compressibility and possible 

failure of reservoir rock 
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Compaction drive 
Compressibility chart 

Effective rock compressibility xIO 6 
[Change in pore vclume/u nit pore volume/PSI 



O Limestone 
• Sandstone 


Rock compressibility 
(pore volume 
compressibility) 
increases when 
porosity decreases 
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Compaction drive 
Ekofisk 



► Ekofisk: a giant offshore field in the Norwegian North Sea 

• Fractured chalk section at about 10,000 ft SS with pay height of about 1,000 ft at crest 

. STOOIP = 6000 MMstb 


► Ekofisk and compaction drive 

• Reservoir was overpressured by about 2,000 psi and initial chalk porosity was very high: 
25 to 48% => pore compressibility increased from 6.10' 6 psi' 1 to a maximum of 10' 4 psi' 1 
providing 30% of total drive energy 

► Ekofisk and subsidence 

• First well drilled in 1969 

• Full field production began by 1974 

• By 1984, due to compaction, seabed subsidence mounted up to 10 feet => platform legs 
had to be cut and new legs section were added up to 20 ft !!! 


Drive Mechanisms - Primary Recovery 


!FP Training 


Compaction drive 



Ekofisk - 2 


GEV: Chalk fields - Impact of subsidence 


t Surface facilities to be repfaced 

• Whiskey WHP replaced by SS Wat. inf. 2010 

• vaihaii VRD 2010/201 1 (Living Quarters and 
Process) 

• Ekofisk Living Quarters 2013 

- Eldtisk II 2015 (WHP + LQ+ Process) 

» Wells to be constantly redrilled (average 
well life 10-15 years) 


► 4D effect visible on seismic: LoFS 



it: 






© 


Ekofisk II replaced 
in 1998 the old 
Ekofisk I facilities 


Total 
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Solution gas drive 



General principle 

► In this type of reservoir, the main 
source of energy results of free 
solution gas liberating from the oil 
and subsequent expansion of this 
solution gas as the pressure is 
further reduced. 

► Oil displacement by the expansion 
of the gas liberated from solution as 
pressure decrease below the bubble 
point pressure 

• As pressure further decreases, more 
gas is released from solution 

• Fluids and rock compressibility effects 
can be neglected vs. expansion of the 
liberated gas => pore volume is 
typically considered to be constant 
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Solution gas drive 
Performances 



► Typical reservoir performances 

• Solution gas drive reservoirs exhibit typically rapid 
reservoir pressure decline and, correspondingly, 
rapid oil production decline but less than in the 
case of natural depletion; after blow-down, the 
decline is very fast for both pressure and 
production 

• GOR quickly reaches a maximum value, before 
declining rapidly (reservoir blow-down) 

• Oil production gets more and more difficult as 
more and more solution gas is liberated 

• Indeed, as a result of the increase in gas saturation in 
the reservoir, there is an unfavorable evolution of 
relative permeabilities and oil flow rate decreases 
while gas flowrate increases 

• Furthermore, the oil viscosity increases since the 
lighter components of the oil get out 

• The recovery factor typically ranges from 5 up to 
25% 
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Solution gas drive 
Performances - 2 



► Typical reservoir performances for undersaturated oil expansion followed by 
solution gas drive 

• The change in the slope of the pressure curve indicates the bubble point 

• GOR is constant as long as P > P b 

• As soon as P < P b solution gas begins to be liberated, but it is not mobile as long as 
S g < S gc where S gc is the critical gas saturation 

• S gc is typically 5% but can go up to 10% for vuggy carbonates 


GOR is slightly decreasing 
since liberated gas is not 
mobile yet 




Sg > Sg C , liberated gas is mobile and GOR 
increases 
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Gas Cap Drive 
General Principle 

► The pressure decline associated with the 
production of oil allows the gas cap to 
expand thus providing further energy for 
additional oil production 

► Gas cap drive effectiveness depends on 

• A large (initial or secondary) gas cap 

• The initial pressure in the oil column @ GOC is 
equal or very close to P b and solution gas drive 
occurs together with gas cap drive 

• The liberated solution gas may form a 
secondary gas cap provided we have high 
vertical permeabilities and a good reservoir 
homogeneity 

• Good contact between oil pool and gas cap 

• Good gravity segregation characteristics 

• Thick oil zone - High dip angle 

• High permeability - especially vertical 
permeability 
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Gas cap drive 
Performances 




► Typical reservoir performances 

• A gas cap drive exhibits typically a slow reservoir 
pressure decline and a slow oil production decline 

• GOR rises slowly and progressively but there is no 
maximum (no reservoir blow-down) 

• Producing slowly will allow the liberated solution 
gas to rise up to the gas cap (and possibly form a 
secondary gas cap) and will help prevent coning 

• Coning is due to the higher mobility of gas (and water) 
compared to oil 

• Gas coning is very difficult to avoid and almost 
impossible to control => we need to perforate lower 
(use Neutron Behind Casing log to know where is the 
GOC) or finally shut the well down 

• The recovery factor is generally higher than for 
solution gas drive reservoirs, 20% up to 40%, and it 
highly depends on the vertical permeabilities 

• Take care not to produce the gas cap! 
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Water drive 
General principle 

► Energy is provided by a water influx from a 
nearby aquifer into the reservoir resulting 
in pressure support and sweeping 

• The pressure decrease in the reservoir under 
production is transmitted to the water zone 

• Hence, the water zone shows a pore volume 
reduction and a water expansion thus leading 
to water production into the oil zone 

► The water drive effectiveness depends on 
the aquifer properties rather than the 
reservoir properties 

• In the short times, the aquifer transmissibility 
that rules the time needed for the aquifer to 
react (i.e. produce water into the reservoir) 

• In the long times, the aquifer size that rules 
the volume of water produced in the reservoir 
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Water drive 

Types of aquifers - infinite aquifers 

► Degree of pressure maintenance 

• Active/ partial / limited water drive 

► Active water drive relates to very active and 
effective aquifers with water encroachment 

• Typically swallow aquifers connected to surface or 
aquifers connected to mountain basement => infinite 
aquifers 

• Infinite aquifers do not exist actually and they refer to 
formations where pressure changes at the water/oil 
contact never reach the outer boundary => outer 
boundary pressure remains constant and equal to the 
initial aquifer/reservoir pressure 

• Bounded aquifers refer to a formation where pressure 
changes at the water/oil contact reach the outer boundary 
i.e. the aquifer outer limit is affected by the water influx 
into the oil pool => outer boundary pressure changes 
with time 
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OIL FIELD ASSOCIATED WITH AN INFINITE AOUIFER 

Rain 








e. MAPVEW 
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Water drive 

Types of aquifer geometries 

► Radial bottom-water drive 

• Radial flow from a large contact area below the 
reservoir 

• Significant vertical flow: driving parameter is k v 


► Radial edge-water drive 

• Radial flow from the side of the reservoir 

• No significant vertical flow : driving parameter is k 


► Linear edge-water drive 

• Linear flow from one side of the reservoir with a 
constant cross-section area 
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Water Drive 
Performances 



► Typical reservoir performances 

• The total fluid rate is generally supposed 
to be constant i.e. water influx may be 
equal to the total reservoir production 
rate thus providing a full pressure 
support (reservoir pressure is 
maintained) 

• If reservoir pressure is kept above 
bubble-point, GOR remains constant 

• A water drive reservoir generally shows a 
steady increase in Water-Cut (and WOR) 
and a decrease in oil production 

• The recovery factor is generally high and 
itcan reach value as high as 40 to 60 % 
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Combined drive mechanisms 
General principle and gravity drainage 



► Several drive mechanisms may exist 

at the same time in the reservoir 

• Depletion drive and a weak water drive 

• Solution gas drive with a small gas cap 
and a weak water drive... 

• The relative importance of each 
mechanism may vary along the life of the 
reservoir 


► Case of gravity drainage 

• Generally not the main drive mechanism 
but rather an additional drive mechanism 

• Gravity drainage effectiveness is related 
to fluid segregation in the reservoir 

• High contrast of gravity between the fluids 

• High reservoir dippage 

• High vertical permeabilities 

• Low production flowrates 

• It may result in very high Rf, up to 60% 
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Primary recovery of oil reservoirs 


0 ^ 'N 

► Primary recovery refers to natural drive mechanisms that can be split in 
two groups: 

• Depletion drive: the reservoir is not in contact with a large active aquifer 

• Monophasic fluid expansion & pore volume compaction 

• Solution gas drive 

• Gas cap drive 

• Water drive: the reservoir is in contact with a large active aquifer that maintains 
pressure 



► From the study of the drive mechanisms, we can get material balance 
equations that are a basic tool for reservoir engineers to analyze the 
reservoir behavior and performance 
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Undersaturated oil expansion and compaction drive 



\ 

► The natural depletion of undersaturated oil reservoirs is the simplest drive 
mechanism 

► It is based on expansion of volumes of oil and connate water and 
compaction of pore volume due to the pressure drop in the reservoir 

► Due to the low compressibilities of fluids and pore volume, the expected 
Rf is low, typically a few % 

► Alternatively, in some cases (overpressure and/or high porosity), the 
compressibility of pore volume can reach high values causing reservoir 

compaction and leading to higher values of Rf, typically 20% 

► However, compaction drive may cause surface subsidence 
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Solution gas drive 



0^ 

► As the pressure drops below the bubble point pressure, solution gas is 
liberated from the oil in the reservoir, it expands and helps produce the oil 

► Solution gas drive reservoirs typically show rapid reservoir pressure decline 
and, correspondingly, rapid oil production decline albeit less rapid than for 
undersaturated oil expansion drive. 

► GOR rises rapidly from initial value to a maximum value before declining 
rapidly (reservoir blow-down) 

► The recovery factor typically ranges from 5 to up to 25 % 
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► The pressure decline associated with the production of the oil causes the 
gas cap to expand 


► To be effective a large gas cap is necessary (initial or secondary gas cap 
formed from solution gas) 

► An efficient gas cap drive mechanism exhibits typically slow reservoir 
pressure decline and a slow oil production decline while GOR rises slowly 
and progressively (no reservoir blow-down) 

► Rf is generally higher than for solution gas drive reservoirs, typically 20- 
40%, and depends highly upon the vertical permeabilities 

► One shall ensure not to produce the gas cap and to limit gas coning at the 
producers 


v I 
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Water drive 



0^ 

► The pressure primary source of energy is supplied by water influx (from an 
adjacent aquifer) into the reservoir 

► The energy for this water influx comes from the compaction of the pore 
volume and the expansion of water in the aquifer 

► Water drive effectiveness depends on the aquifer properties and not on 
the reservoir ones; the two key parameters are the aquifer transmissibility 
in the short times and the aquifer size in the long times and 

► The total fluid rate remains generally constant, with a sharp increase in 
Water-Cut while GOR remains constant as long as pressure is maintained 
above P b 

► The recovery factor typically ranges from 30 to 60 % 
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Principles of the 
material balance 


Principles of the material balance 
Calculation methods for production/reserves forecast 

► Simplified methods 

• Material balance: single cell model with simplified expression of the recovery factor 

• Decline rate analysis: history match of flowrate (wells, field) in order to forecast reserves 

► Advanced methods 

• Numerical models: building a model in order to mimic the behavior of the actual 
reservoir 

• Numerical models allow to handle heterogeneities and related uncertainties at all scales 

• Numerical models allow to handle complex fluid behavior 
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Introduction to primary recovery 
Material balance vs. numerical simulation 



The geological/reservoir models are not 

required 

• Production data 

• Pressure data 

• PVT data 



Accumulation 

Identification of production mechanisms 
Production forecast (cumulative) 


The geological/reservoir models are 

needed 

• Production data 

• Petrophysical properties 

• PVT data 

Q 

Volume of hydrocarbons in place 
Production mechanisms 



Pressure value forecast 
Saturation value forecast 
Production forecast (flowrate and cumulative) 
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Principles of the material balance 
General principle 



► Material Balance = balance regarding the volumes of fluids in the reservoir 

• Assuming the reservoir pore volume is 100% filled with fluids, in reservoir conditions: 

• Initial Vol. fluids = Remaining Vol. fluids + Change in pore volume + Vol. of injected or entered 
fluids 

Or 

• Net withdrawal fluids volume (withdrawal - injection) = Expanded vol. of fluids in the system + 
Vol. of cumulative water influx 






A kind of "equation of continuity" for the reservoir/part of the reservoir over a finite 
time interval => the simplest reservoir simulation model (1 cell!) 

It requires a well established pressure regime over the reservoir or part of the 
reservoir 



After oil production has begun 



Producer 


pore compaction and connate water expansion 

water intliiZI 


Pj/ B oi , B gi , R si 
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Principles of the material balance 
Conditions for applying the material balance 

► There are no 'sufficient' conditions for the meaningful application of the 
material balance to a given reservoir, but there are two 'necessary' 
conditions that must be satisfied 

• There should be adequate data collection (production data/pressure data/PVT), both 
in quantity and quality, otherwise the attempted application of the technique becomes 
meaningless. 

• Production data are typically average parameters e.g. average GOR (not instantaneous GOR) => 
they must be calculated 

• Pressure data classically refers to reservoir static pressure not bottom hole flowing pressure => 
it must be measured or evaluated 

• It must be possible to define an average pressure decline trend for the reservoir or 
part of the reservoir 

• The pressure dependent PVT parameters and Ap are evaluated as a function of time using this 
decline trend 

• The pressure decline trend may be defined separately on various compartments of the 
reservoir => The material balance may be applied separately to the various compartments 
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Principles of the material balance 
Accumulation, production and injection 

► Evaluation of accumulation (to be compared with volumetric methods) 


Accumulation 

N 

G 

w 

Cumulative 

production 

N p 

G p 

W P 

Cumulative 

injection 


G lnj 

W lnj 


All volumes in standard conditions 


► We need to know the Formation Volume Factors to convert from surface to 
reservoir conditions 
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Gas Water 


Oil 






Principles of the material balance 
Fluid volumes before production 




Surface 

Volumes 

(standard m 3 , 
stb Standard Barrels 
or scf Standard Cubic Feet) 


Subsurface Volumes 

(reservoir barrels 
reservoir m 3 ) 
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Principles of the material balance 
Fluid volumes after production 



SCF 


STB 


STB 



APV 


Surface volumes 
(standard m 3 , stb or scf) 


Compaction 


Water entries (externally derived water) 


Subsurface volumes 
(reservoir m 3 or barrels) 


Released Solution Gas 
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Symbols and units 



Symbols UNITS 




US units 

Metric 

Mg 

Gas viscosity 

cP 

cP 

Mo 

Oil viscosity 

cP 

cP 

Bg 

Gas formation volume factor at pressure P 

Bbl/Mscf 

vol/vol 

Bgi 

Initial gas formation volume factor 

Bbl/Mscf 

vol/vol 

B 0 

Oil formation volume factor at pressure P 

Bbl/stb 

vol/vol 

B oi 

Initial oil formation volume factor 

Bbl/stb 

vol/vol 

B, 

Total hydrocarbon formation volume factor 

Bbl/stb 

vol/vol 

B,i 

Initial total hydrocarbon formation volume factor 

Bbl/stb 

vol/vol 
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Symbols and Units - 2 



Co 

Oil compressibility 

psi' 1 

bar' 1 

c p 

Pore compressibility (c f ) 

psi' 1 

bar' 1 

c w 

Water compressibility 

psi' 1 

bar' 1 

p 

Pressure 

psi 

bar 

Pb 

Bubble point pressure 

psi 

bar 

Pi 

Initial pressure 

psi 

bar 

Rl 

Liberated gas 

scf/bbl 

m 3 /m 

Rs 

Solution gas at pressure P 

scf/bbl 

m 3 /m 

Rsi 

Initial gas in solution at Pj 

scf/bbl 

m 3 /m 


lbbl = 5.615 cuft 
lbbl = 0.159 m 3 

P std = 14.7 psia T std = 60 °F 

P std = 1.01325 barsa T std = 15 °C 

1 bar = 

14.5 psi 
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Principles of the material balance 


\ 

► The material balance expresses a conservation of the volumes of the fluids 
for the reservoir / part of reservoir for a finite time interval 

• In order to be valid, a well established pressure regime shall exist in the reservoir or 
part of the reservoir under consideration 

► Since the accumulations and produced volumes are expressed in surface 
conditions we need to use fluid properties like the formation volume 
factor in order to convert from surface to reservoir conditions 

► Compressibility effects are of great importance since they are the basis for 
volume conservation: oil expansion, gas expansion, connate water 
expansion and pore volume reduction (compaction) 
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Note 
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Simplified MBE for undersaturated oil expansion 
► Considering that there is no aquifer activity (i.e. no water entry in the 


reservoir), no gas cap and no liberation of solution gas and disregarding the 
compressibility of the connate water and of the pore volume 

► It can be stated that the volume occupied by the oil remains constant: 

volume of expanded remaining oil = initial volume of oil 

(N - N p )B 0 = NB ol 
Rf = N P /N = 1 - BJB 0 

► However water and pore volume compressibilities have to be taken into 
account to better estimate the recovery factor 

• Water and pore volume compressibilities are generally of the same order of magnitude 
than oil compressibility 

IFP Training I 57 

Drive Mechanisms - Primary Recovery — ' 


Material balance equations for oil reservoirs 
Detailed MBE for undersaturated oil expansion 



► Still considering that there is no water entry, no gas cap, no liberation of 
solution gas 

Produced volume = Increase of volume of oil (expansion) 

+ Increase of volume of water (expansion) 

+ Decrease of pore volume (compaction) 

► For a pressure drop A P from P t to P with P>P b (undersaturated oil): 


• 

Produced oil volume: 

Np-Bo 


• 

Increase in volume of oil: 

{y p ■ s 0 ■ c o) ■ 

A P 

• 

Increase in volume of connate water: ( V p ■ S wi ■ c w ) 

■A P 

• Decrease in pore volume: 

where V p is the pore volume 

(Vp ■ Cp) ■ 



N P B 0 

= V p AP(c 0 S 0 + c w S wi 

+ C p ) 


Drive Mechanisms - Primary Recovery 


!PP Training 


58 


Material balance equations for oil reservoirs 
Detailed MBE for undersaturated oil expansion and Rf 



► Introducing equivalent compressibility: 

C e ( "I” Cw^wi "I” Cp)/ S q Cf/S 


o 


► The material balance equation (MBE) writes: 

N p B 0 = c e S 0 V v AP 


with S Q V p = NB oi 


N V B 0 = NB oi c e AP 


Rf = Np/N = 




c e AP 


► From this equation, compressibilities of oil, water and pore volume being 
generally pretty low, we can infer that the calculated recovery factor will be 
low, typically a few %. 
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Material balance equations for oil reservoirs 
Relationship between c Q and B 0 above bubble point 



► By definition, c 0 is the fractional change in oil volume with a unit change in 
pressure. 

• This can also be expressed as a function of the variation in B 0 above bubble point 

-1 dV 

C °~~V~'dP 

hence 

-1 dB 0 
C ° ~ ~B~ 0 ~dP 
or 

c„ = -KB 0 - B ot )/B oi ]/(P - PO 


► Hence, above the saturation pressure, c 0 can be derived from the slope of the 
line B 0 (P) 
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Material balance equations for oil reservoirs 
Using the undersaturated oil expansion MBE 

► First use: 

• Knowing N, A P and the fluid properties, estimate N p and Rf 

• In general, Rf < 5% leading to a secondary recovery scheme (e.g. water injection) 

► Second use: 

• Knowing N p and A p (from production history), calculate the accumulation N and 
compare to the previous estimate from volumetric studies 

• If N is found to increase over time, it means that actual pressure drop in the reservoir is 
lower than the measured one => we may suspect a pressure support from an active 
aquifer and the drive mechanism may be natural water drive 

• If N is found much lower than the estimate coming from volumetries, it means that only 
one part of the reservoir is drained 

• Compartmentalization => we may think to an additional development (new wells) 

• Heterogeneities in permeability => possible reconnection of the missing accumulation (slow) 
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Material balance equations for oil reservoirs 
MBE for solution gas drive 

► Considering there is no water entry and no gas cap 

► Solution gas drive 

• Reservoir pressure drops below bubble point pressure and free solution gas is liberated 
in the reservoir 

• The main mechanism is expansion of the released solution gas 

• Therefore the expansion of connate water can be neglected because of the respective 
value of connate water compressibility and gas compressibility 

• Accordingly, the compaction of the pore volume can also be neglected => pore volume is 
considered to be constant 

• Therefore 

oil volume + liberated gas volume = constant 

• At initial pressure P t (= Pb): V oi , V w at current pressure P : V or , V w> V gf 
Hence 

V()i ~ V OT + V gf 

• Finally introducing R s the solution GOR and B g the gas FVF, the final MBE is given by 

NB ol = (N — N p )B„ + [NR si - (N - N p )R s - C p ]B g 
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Material balance equations for oil reservoirs 
PVT data reminder 



> Reminder of variation of main reservoir oil properties with pressure 
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Material balance equations for oil reservoirs 
MBE for solution gas drive and Rf 

► Introducing a new term in the previous equation: average GOR noted R p 

Rp = Gp/Np 
Gp — RpNp 

► The MBE for solution gas drives can be expressed as: 

N p [B 0 + (R p - R s )B g ] = N[(B„ - B ot ) + (R sl - R S )B„] 

oil production free gas production oil expansion solution gas expansion 


► Therefore the recovery factor can be expressed as: 

Rf = N p /N = [<B„ - B ot ) + (R sl - R s )B g y[B 0 + (R p - R S )B„] 

Hence Rf = f(l/Rp) : the recovery factor is a function of the inverse of the average GOR 


► MBE shows that as much gas as possible must be kept into the reservoir 
• Production flowrate should be low 
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Material balance equations for oil reservoirs 
Total formation volume factor and new MBE 



> Reminder about the variation of the FVF (Formation Volume Factor - B) with 
pressure 



> Introducing the total FVF in the solution gas drive MBE 

N P [B 0 + ( R p - R s )B g ] = N(B t - B ti ) 

Rf (B( Bfi)/[B 0 + (Rp 
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Material balance equations for oil reservoirs 
MBE for gas cap drive 

► Still considering there is no water entry and neglecting the connate water and 
pore volume compressibilities 

► Gas cap drive 

• The main mechanism is the expansion of the gas cap that comes in addition to the 
expansion of the liberated solution gas 

• Hence, we add a term expressing the gas cap expansion to the solution gas drive MBE 

• Introducing the parameter m = Volume of gas zone / Volume of oil zone: 

m = GB gi /NB oi 

• Therefore 

G = mNB oi /B gi @Pj and GB g = mNB oi B g /B gi @P « Pj 

• Hence the gas cap expansion is given by: 

GBg - GB gi = mNB oi B g /B gi -mNB oi = mNB oi (B g /B gi - l) 

• Therefore the final MBE for the gas cap drive is: 

N p [B 0 + ( R p — R s )B g ] = N[(B a — B oi ) + (R si — R s )B g + mB oi (B g /B gi — l)] 

oil production free gas production oil expansion solution gas expansion gas cap expansion 
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Material balance equations for oil reservoirs 
MBE for gas cap drive - 2 



> The previous equation is quite complex and does not provide a clear picture 
of the principles involved in the gas-cap drive mechanism 

> However, because of the gas-cap expansion, the pressure decline is less 
severe than for solution gas drive and oil recovery is greater 

> If not neglecting water and pore compressibility, the final equation is: 


Np\B 0 + i^Rp Bs^Bg] 


= N 


(Bo ~ B 0 i) + ( R S i ~ R s )B g + mB oi (B g /B gi — l) 


connate water expansion and pore volume compaction 


within oil zone AND gas cap 
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Material balance equations for oil reservoirs 
MBE for natural water drive 



► Natural water drive implies water entry from the aquifer into the reservoir 

► The corresponding volume of water has to be added to the corresponding 
equation without water drive 

► Assuming undersaturated oil expansion (for sake of simplicity) and 
considering all fluids and pore volume compressibilities 

Oil production = Oil expansion + Connate water expansion + Pore volume compaction 
+ Aquifer expansion => water entry W e 
- Water production => water production^ 

Therefore 

N p B 0 = NB oi c e AP + W e - W p B w 


Assuming no water production, the recovery factor is: 


B 


Rf = N v /N = ~^c e AP + 


w, 

NB, 


' o 

In any case, we need a water influx model in order to calculate W e 
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Material balance equations for oil reservoirs 
Water influx models 

► Water influx models 

• Pot aquifer model 

• Schilthuis steady-state model 

• Hurst modified steady-state model 

• Hurst and Van Everdingen unsteady state model 

• Radial edge-water drive 

• Radial bottom-water drive 

• Linear edge-water drive 

• Carter-Tracy unsteady-state model 

• Fetkovitch method 

• Radial aquifer 

• Linear aquifer 

► Water influx calculations are generally very uncertain since aquifer 
characteristics are poorly known and heavily rely on historical production 
data in order to indirectly estimate the aquifer characteristics 



Drive Mechanisms - Primary Recovery 


!FP Training 


69 


Water influx models 
Pot aquifer model 



► Pot aquifer model = the simplest water influx model 

• The pressure drop due to the production in the reservoir causes the 
aquifer to expand and the water to flow into the reservoir 

• Writing the corresponding MBE, we get: 

W e = W,c,AP = W,(c„ + Cf)-f- (P, - P) 
where W t is the initial volume of water in the aquifer 
c t is the aquifer total compressibility 
Pi is the initial aquifer/reservoir pressure 
P is the current reservoir pressure (@WOC or @GWC) 

/ is the partial encroachment angle 

• This model is valid for small aquifers only (i.e. r e « r 0 ) where the 
pressure drop at the aquifer/reservoir boundary transmits almost 
instantaneously through the entire aquifer 

• In the case of large aquifers, a time dependent equation is needed: 
the aquifer needs some time to react to the pressure change in the 
reservoir and to begin to produce inside the oil pool 
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Water influx models 

Hurst & Van Everdingen unsteady-state model 



► The equation of water flowing from the aquifer into the reservoir is derived 
from the equation of oil flowing from reservoir into the well 

• In the latter case, the solution in pressure writes: 

, N 2nkh , . 

P D< t D)=—(P- P wf) 

where P D (t D ) is the dimensionless pressure function corresponding to the CTR case 


• In the case of determining water influx, we are much interested in calculating the influx 
flow rate 

• Hurst and Van Everdingen solved the radial diffusivity equation for the reservoir-aquifer 
system and derived the constant terminal pressure (CTP) solution of form: 

q ° ( ^ 2nkhAP 

where qo(.t D ) is the dimensionless influx rate evaluated at r D = 1 

• Vd^d) describes the change in rate from 0 to q due to a pressure drop A P applied at 
the outer reservoir boundary r Q at time t = 0 
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Water influx models 

HVE model - Solution for edge-water drive 



► The solution is generally expressed in terms of cumulated influx by 
integrating previous equation 

► Therefore HVE solution writes: 

foQdt = f* D q D (t D )-^dt D ■' * k 


p 


2nkhAP 


with t D = 


&\ic t r 0 


pW e 

2nkhAP 


= W eD (t D ) 


<P\j£tTo^ 

k 


where W eD — f Q D qD(t D )dt D is the adimensional water entry 

► In general 

W e = 2nf0hc t r o 2 ■ A P ■ W eD (t D ,r eD ) = U ■ AP ■ W eD (t D ,r eD ) 

where 

U = 2nf&hc t r 0 2 is the aquifer constant 
WeD is determined from tables as a function of t D and r eD 
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Water influx models 

HVE model - Solution for radial edge-water drive 



• Aquifer constant 
U = 2nf0hc t r o 2 

U = 1.119 f0hc t r o 2 

• Dimensionless time 

_ fc t 
D 0\ic t r a 2 


Darcy units (cc/atm) 
Field units (bbl/psi) 

Darcy units (t-seconds) 


t D = c te — — Field units 

0\ac t r 0 



0 ° 

360 " 


— Radians 
2k 


With c te = 0. 000264 t-hours 
c te = 0. 006328 t-days 
c te = 2. 309 t-years 
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Water influx models 
HVE model - Practical use 



► The former equation gives the expression of the water influx due to an 
instantaneous pressure drop at the reservoir/aquifer boundary 


► In practice, we get a continuous pressure decline at the reservoir/aquifer 

boundary 

• In this case, the continuous pressure decline is approximated by a series of discrete 
pressure steps A 

• Then, for each pressure drop step, the corresponding water influx can be calculated by 
using the previous equation 

• The final cumulative water influx is given by superposition of each water influx, by 
respect to time 


W e (P l .i) = u'S' n W eD (t D ,r cD )APi 

Z— i/=l 
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Water influx models 

HVE model - Practical solving method 



► Solving HVE model procedure 

• Calculate r eD = r e /r Q and U 

• Discretize the pressure history 

• For each time step 

• For each pressure step 

• Calculate the corresponding value of adimensional 
time t D 

• Calculate the corresponding value of W eD from the 
tables and/or curves 

• Sum up the contribution of each pressure drop and 
find the corresponding value of water entry W e 

=> HVE method is an additive superposition process 
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Water influx models 

HVE model - Solution for radial bottom-water drive 

► The HVE solution to the radial diffusivity equation is considered as the best 
aquifer model to date but it is not adapted to bottom-water drive since it 
cannot model vertical encroachments 

• A model was derived by Coats (1962) taking into account vertical flows through a new 
parameter, the permeability ratio: F k = k v /k h 

• A general solution was derived by Allard & Chen (1988) very much similar to the one of 
Hurst & Van Everdingen: 

W e = 2n0hc t r o 2 ■ A P ■ W eD (t D ,z D ) = U ■ AP ■ W eD (t D ,r eD ,z D ) 

where z n = — is the dimensionless vertical distance 
D r 0 jr k 

• The final solution is of the same nature: 

W e (Pj, t) = l/^ W e 0 (t D , r eD , Zd)APi 

• The values of adimensional water entry W eD (t D ,r eD ,z D ) are tabulated but they differ 
from the ones used for HVE solution for the edge-water drive 
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Water influx models 
Carter-Tracy unsteady-state model 



► HVE model provides the exact solution to the radial diffusivity equation but 
the use of superposition principles involves quite tedious calculations 


► Carter and Tracy proposed a technique (1960) allowing a direct calculation of 
the water influx 


• The water flow rate is assumed constant over a finite time interval corresponding to a 
given pressure drop => the cumulative influx at time n can be calculated from the 
previous value at time n — 1 using the equation: 


/ UAp(t D j ) — W e (t z , y -_ 1 ) ■ P' 


"•(**) = + 1 ;X)-;^.(Cn • (t * - 1 °- ) 


Where U is the HVE water influx constant 

Pd(Id) ' s the dimensionless solution of diffusivity equation with constant rate 
boundary condition 

= dP D (.t D )/dt D is its time derivative 
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Water influx models 
Solving Carter-Tracy model 

► Several possibilities: 

• The Pd^d) have been tabulated by Hurst and van Everdingen and can be used to 
directly compute W e (t D ) iteratively 

• If t D > 100 the following expressions can be used: 

P D (t D ) = 0. 5 (In t D + 0. 80907) 

P D ' (t D ) = 1/(2 ■ t D ) 

• It is also possible to use the Fanchi's regression equation: 

P d^d) = a o + a \tD + ^2 In to + #3 (In t/)) 2 

The values of the regression coefficients are given for different values of the ratio r eD = r e / r 0 
(from Dake-The Practice of Reservoir Engineering) 


d 

ao 


£12 

ay 

u 

0.10371 

1.66657 

-0.04579 

-0.01023 

2.0 

0.30210 

0.68178 

-0.01599 

-0.01356 

3.0 

0.51243 

0.29317 

0.01534 

-0.06732 

4.0 

0.63656 

0.16101 

0.15812 

-0.09104 

5.0 

0.65106 

0.10414 

0.30953 

-0.11258 

6.0 

0.63367 

0.06940 

0.41750 

-0.11137 

8.0 

0.40132 

0.04104 

0.69592 

-0.14350 

10.0 

0.14386 

0.02649 

0.S9646 

-0.15502 

00 

0.82092 

-3.68 x 10" 4 

0.28908 

0.02882 
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Material balance equations for oil reservoirs 


0 ^ 

► From the material balance principle, equations can be written for each of 
the main natural drive mechanism 

• These equations confirm the low recovery factor expected for the undersaturated oil 
expansion and the higher values for solution gas drive, gas cap drive and water drive 

• In the case of solution gas drive and gas cap drive, the MBE clearly shows that the gas 
must be kept as much as possible within the reservoir in order to get a high recovery 

• The natural water drive MBE needs water influx in order to predict recovery 

• Several models exist, the better being Hurst and Van Everdingen transient model 

► Using MBE can be two-folds 

• It makes production forecast by calculating Rf between an initial and a final pressure 

• It uses production data to test hypotheses, typically the existence of an active aquifer 

• It uses production data to match the reservoir description parameters 

• Typically the aquifer parameters (especially the size), the size of the gas cap (through 
parameter m) and finally accumulation 

• Caution: MBE matches only connected accumulation 
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Generalized material balance for oil reservoirs 
Most general MBE for oil reservoirs 

► Considering the following general case: 

• A reservoir with an oil pool overlain by a gas-cap and underlain by an aquifer. 

• The oil pool is produced and the production at the surface will consist of oil, gas and 
water. 

• Initially, there is no injection of water or gas. 


► The most general MBE writes 

Np\Bo + (^P — 


= N 


iPo — B oi ) + ( R S i — R s )B g + mBoi^Bg/Bgi — l) 


+ (1 + m)B oi I Cm,5ivC - + Cp MP 


1-5 


WC 


+ w e - W p B w 
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Generalized material balance for oil reservoirs 
Drive index 



► Introducing 

D = N p [B 0 + (fl p — R s )B g ] = N p [B t + (R p — R S i)B g ] 

and dividing by D we get: 

N[(B 0 B 0 i) + — Rg)Bg ] mNB 0 i ^B g /B g i l] 

= D + D 


(1 + Wl)NB 0 i f c w S wc + Cp\ W e W p B w 

+ ^ — — - AP + — ^ — 

D y 1 — S wc ) D 


1 = DDI 

Depletion 
Drive Index 


+ 


SDI 

Segregation 
Drive Index 


+ CDI 

Formation / Connate Water 
Compressibility Drive Index 


+ WDI 

Water 
Drive Index 
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Generalized material balance for oil reservoirs 
Interpretation of the drive index 

> In the previous equation: 

• The numerators of the fractions represent: 

• the expansion of the initial oil zone (included released solution gas), 

• the expansion of the gas cap, 

• the compaction of the pore rock volume / expansion of connate water 

• the net water influx. 

• The denominator represents the volume of hydrocarbons produced, expressed at 
current reservoir pressure conditions 

• Each of the parameters DDI, SDI, CDI and WDI, characterizes one of the production 
mechanisms of the reservoir 

• Undersaturated oil expansion and solution gas drive 

• Compaction gas drive 

• Connate water expansion and compaction drive 

• Natural water drive 

• The evolution with time of the drive index will help identify the prominent drive 
mechanisms and make decisions for the future development 
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Generalized material balance for oil reservoirs 
Injection drive index 

► In the case where there is a fluid injected within the reservoir 

• Secondary recovery (see below) 

• Either water or gas or water and gas 

► In this case, a new additional drive index is added to MBE: Injection drive 
index 

• Located in the right hand side of the previous equation 

G- B ■ • 4- W- B ■ ■ 

IDj ^inj 1 * ginj ' vv inj^winj 

~ D 

• Note that very often B winj - = B w but B gin j =£ B g since the properties of the injected gas 
(after treatment) are different from the properties of the produced gas 
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Generalized material balance for oil reservoirs 
Another expression of the generalized MBE 

► Another way to express the material balance is to state that 

net total fluid withdrawal = expansion of the HC fluids in the reservoir 

+ cumulative water influx 

• The net total volume of fluid withdrawal is given by: 

F = N p {B 0 — R S B g ) + Bg(^Gp — G in j) + B w (W p — 

• The expansion of oil and its original dissolved gas is: 

NE„ = N[(B 0 - B oi ) + (R si - R s )B g \ 

• The expansion of the gas cap is 

NmE g = NmB ot (B g /B gi - l) 

• The pore volume compaction / connate water expansion is: 

( ^ wc^w "f" 

NE fiW = N(1 + m)B oi f 15 M A P 
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Generalized material balance for oil reservoirs 
Another expression of the generalized MBE - 2 



► The new form generalized material balance equation writes: 


F — N(E 0 + mEg + Efw) + W e 


F = NE t + W e 


Where: 


Ef — E 0 + mE g + E f w 


Eo — (B a — Boi ) + (E si — R s )B g 
Eg — B 0 i{Bg/Bgi — l) 

E fjW = (l + m)B ot (^^)AP 


is the total expansion factor 
is the oil expansion factor 
is the gas expansion factor 

is the connate water expansion / PV compaction factor 
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Generalized material balance for oil reservoirs 
Material balance analysis - Havlena-Odeh method 



► In order to analyze the generalized MBE, Havlena and Odeh proposed a 
methodology that used plotting techniques to try to show a linear behavior 
from the equation 

• Havlena D. & Odeh A. S., The Material Balance as an Equation of a Straight Line, Journal 
of Petroleum Technology, Aug. 1963 


► In the general case by plotting F as a function of E 0 + mE g + Ef w we may 
have a straight line of slope N and intercept W e 

► However, interpreting this plot can be quite difficult especially if N, m and 
the aquifer model are uncertain 


► In most of the specific cases, the generalized MBE simplifies and may lead to 
simple plots 
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Generalized material balance for oil reservoirs 
Havlena-Odeh method - Solution gas drive 



► Assuming 

• No gas cap and no aquifer 

• Pressure decreasing below bubble point pressure 

• Connate water and pore volume/rock compressibilities neglected 


► The generalized MBE writes: 

• F = NE 0 => F/E 0 = N 

• => by plotting F = f(E 0 ) we may have a straight line of slope N 



Generalized material balance for oil reservoirs 
Havlena-Odeh method - Gas cap drive (1) 

► Assuming 

• No aquifer 

• Pressure decreasing below bubble point pressure 

• Connate water and pore volume/rock compressibilities are neglected 

► The generalized MBE writes: 

• F — N(E 0 + mE g ) 

• => by plotting F = f{E 0 +mE g ) we may have a straight line of slope N 

This plot can be used to check different hypotheses on the 
value of m 

The correct value of m is obtained when the plot shows a 
straight line of slope N passing by the origin (assuming N is 
known) 
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Generalized material balance for oil reservoirs 
Havlena-Odeh method - Gas cap drive (2) 



► Same assumptions 

► The generalized MBE writes: 

• F = N(E 0 + mE g ) => F/E 0 = N + mN^- 

• => by plotting F/E 0 = f{E g /E 0 ) we may have a straight line of slope mN and intercept 
N 
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This plot can be used to check hypothess on m 
and N 

When the plot shows a straight line, the 
intercept with the Y axis gives N and the slope 
gives mN 
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Generalized material balance for oil reservoirs 
Havlena-Odeh method - Water drive 



► The generalized MBE writes: 

F = NE t + W e => F/E t = N + W e /E t with E t = E 0 + mE g + Ef W 

=> by plotting F/E = f(W e /E t ) we may have a straight line of slope 1 and intercept N 


Note that E t expression may simplify depending on the drive mechanisms to be considered 



This plot is very useful to confirm the water influx model 

Indeed, aquifers, even analytically, are characterized by 
various parameters: type, angle of encroachment, external 
radius, petrophysics etc. all uncertain 

The appropriate aquifer model will show a straight line of 
slope 1 and intercept N ; any deviation is due to incorrect 
geometry or incorrect size of the aquifer 
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Generalized material balance for oil reservoirs 
Havlena-Odeh method - Aquifer and gas cap 



► Assuming 

• N and m are known 

• Connate water and rock compressibilities are neglected 

► The generalized MBE rewrites: 

F = N(E 0 + mE g ) + W e => F/(E 0 + mE g ) = N + W e /{E 0 + mE g ) 

=> by plotting F/(E 0 + m Eg) = f(W e /(E 0 + mEg)) we may have a straight line of slope 
1 and intercept N 
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This plot is used in the same way than the previous one 
in order to confirm the water influx model 
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Generalized material balance for oil reservoirs 
General case - Campbell plot 



► In the most general case, MBE writes: 

F — N(E 0 + mE g + Ef W ) + W e 

that is F = NE t + W e => (F - W e )/E t = N 

=> by plotting(F — W e )/E t = f(N p ) we may have a straight line of slope 0 and intercept N 
i.e. of equation y = N 


(F-W e )/E t 


/ 

* 

X 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

1 

strong aquifer 

1 s N 

,/ v 
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moderate aquifer 


? volumetric depletion 

weak aquifer 



N 


To have a good match, all the points should be on the 
line y = N ; deviation from the horizontal line gives an 
indication on the way to modify the aquifer 

The Campbell plot is widely used since it gives good 
feelings about the aquifer characteristics 
In general, the plot is built without any aquifer in the 
model to check which type of aquifer should be added: 

- strong aquifer => constant pressure boundary 

- moderate aquifer => infinite aquifer 

- weak aquifer => no flow boundary 

- no aquifer => volumetric depletion 
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Generalized material balance for oil reservoirs 



\ 

► Oil field generalized material balance equation is quite complex and 
difficult to solve in theory (finding N, m and the aquifer model parameters) 

► Havlena and Odeh methodology proposes a simplified analysis based upon 
plotting and trying to find a linear behavior 

► These plots allow to confirm the value of the original oil in place N; they 
also allow to check the presence of an active aquifer and, if any, to match 
the water influx model 

► In the most general case, Campbell plot is widely used to check the 
presence of an active aquifer and to determine its characteristics 

► The Drive Index is an additional method allowing to analyze the 
generalized MBE and to compare the relative importance of the various 
production mechanisms one to each other, thus providing a guide for 
further development of the field 
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Note 
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Gas reservoirs 
Reminder on gas reservoirs 



► Hydrocarbons are generally classified as either oil reservoir or gas reservoirs, 
depending on the state of the reservoir contents at discovery. 

• For gas reservoirs, the different gas definitions refer to the PVT behavior of the gas, 
which depends on the initial reservoir gas composition, initial reservoir pressure and 
temperature. 

► A gas reservoir is a reservoir in which hydrocarbons remain in the gaseous 
state throughout the life of the reservoir. 

• Dry gas or wet gas reservoirs 

• A reservoir in which part of the gas condensates in-situ is a gas condensate reservoir 



Gas reservoirs 

Gas reservoir drive mechanisms 



► Gas reservoirs are produced mainly through monophasic natural depletion 


► Gas compressibility is much higher than liquids/pore volume compressibilities 

• We expect higher values of recovery factor than for oil reservoirs, typically above 90% 

• The lower the abandon pressure, the higher the recovery factor 

► Main concern of gas reservoirs: aquifer activity 

• In case there is an active aquifer, the reservoir will be flooded by water 

• Despite the difference in mobility in favor of the gas, a significant part of the gas will be 
trapped by the water behind the front 

• We expect a much lower recovery factor than without aquifer activity, typically in the 
range of 60-70% 

• In this case, recommendation is to produce as fast as possible before flooding of the 
reservoir 
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Gas reservoirs 

MBE for gas reservoir without water entry 



► Assuming no water entry and neglecting connate water and pore volume 
compressibilities compared to gas compressibility 

► When pressure drops, the volume occupied by the gas under reservoir 
conditions does not change i.e. 

initial gas accumulation = gas accumulation @(P,T) 


GBgi — (G Gp^Bg 

Gp = G{\ — Bgi/Bg) 


Where 

G is the initial accumulation of gas in standard conditions 
G p is the cumulative gas production in standard conditions 
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Gas reservoirs 

Reminders about equations of state (EOS) 



► Equations of state: f(P, V, T, ri) = 0 

► Ideal gas law (Mariotte, 1650): 

PV = nRT 

► Equation of state for real gas: 

PV = ZnRT 

Where Z is the compressibility factor: a measurement of the difference in the real gas from 
the ideal gas ; when P -> latm then Z 1 

► Cubic equations of state 

• Van derWaals (1873) P = RT/(V - b) - a/V 2 

• Redlich-Kwong (1949) P = RT/(V - b) - a/[V(V - b)T^ 2 ] 

• Soave-Redlich-Kwong (1972) P = RT /(V — b) — a(T)/[V(V + b )] 

• Peng-Robinson (1976) P = RT/(V - b) - a(T)/[V 2 + 2 bV - b 2 ] 

• Coefficients a and b are matched using PVT package typically from results coming from 
CVD experiments (cf. PVT studies) 
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Gas reservoirs 
Expansion factor and FVF 

► Gas expansion factor: 

E = 1/B g 



> B g is the gas FVF i.e. ratio of the volume of n moles of gas at reservoir 
conditions divided by the volume of n moles of gas at standard conditions 

@standard conditions: = Zo n RT 0 = nRT 0 

(©reservoir conditions: PV = ZnRT => B g = Z ■ P Q /P ■ T/T 0 

(©initial conditions: P^i = ZinRTi => B gi = Zi ■ P 0 /Pi m Ti/T 0 

► Assuming that the reservoir temperature remains constant (isothermal 
transformation in the reservoir), we get: 

^ *_gi _ Zj P 
~B g ~~P i Z 
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Gas reservoirs 



Recovery factor for gas reservoir without water entry 
► The recovery factor writes: 



i 


► Recovery is a function of: 

• The original reservoir pressure 

• The current reservoir pressure 

• The gas mixture composition through Z (and Z t ) 

► It is assumed that the reservoir depletes in the same manner everywhere i.e. 
the pressure rate decrease is the same everywhere in the reservoir 

• Corresponding Rf can be very high, up to 90-95% (e.g. Lacq) 
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Gas reservoirs 



Typical behavior of gas reservoirs without aquifer 


► Analysis: plotting P/Z = /(G p ) 



• it is a straight line for a closed reservoir (no water entry) 

• The extrapolation of the line gives the estimate of the accumulation G 


p 


z 



no aquifer W e = 0 


extrapolation gives accumulation 




G 


Produced gas G 
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Gas reservoirs 
Overpressured reservoirs 



► For highly compressible reservoirs (overpressured), the combined effect of 
connate water and pore volume compressibility cannot be ignored at high 
initial pressures and the generalized MBE writes: 

G(1 - c e AP)B gt = (G - G p )B g 

Where c e = 

e 1-Swc 

► The recovery factor writes 

Rf = G p /G = 1 - (1 - c e AP)B gl /B g = 1 - (1 - c e AP) yX^ 

► When the pressure in the reservoir decreases, the effect of c e can be 
neglected again in front of c g 

► In this case, P/Z plot typically displays a change in the slope which can lead 
to very significant overestimate of the accumulation G 

• A change in slope refers to a change in pore volume/rock compressibility when the 
pressure decreases but some overpressured gas reservoir does not show any slope 
change !!! , 
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Gas reservoirs 

Overpressured reservoirs - P/Z plot 

► Analysis: plotting P/Z = /(G p ) 

• Two slopes 

• c g dominates only when pressure decreases enough 
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Gas reservoirs 

MBE for gas reservoirs with water entry 



► Still neglecting connate water and formation compressibilities compared to 
gas compressibility, the generalized material balance equation writes: 

net fluid withdrawal = expansion of the HC fluids in the reservoir 

+ cumulative water influx 


► The net total volume of fluid withdrawal is given by: 

► The expansion of gas is: 


F — G p Bg + W p B w 


GEg - G(Bg - Bgi ) 


g " gi , 

► The generalized MBE for a gas reservoir with water entry is given by: 

T e 


F = GEg + W t 


GpBg + WpB W — G{Bg — Bgi ) + W e 
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Gas reservoirs 

Rf for gas reservoirs with water entry 



► The recovery factor is given by: 

Rf = G p /G = 1 - BgjBg + 1/ GBg ■ ( W e - W p B w ) 


Rf—1 — 


B 


gi 


B 


1- 


W e - W p B w 


9 


GB 


= 1 - 1 - 


W e - W p B w 


gi 


GB 


gi 


Pj/Zj 

P/Z 


► Recovery is a function of: 

• The original reservoir pressure 

• The current reservoir pressure 

• The gas mixture composition 

• Water entries W e 


► The corresponding Rf is lower than without water entry because of the gas 
trapped behind the water front, ranging between 50% and 70% 
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Gas reservoirs 

MBE analysis for gas reservoirs with water entry 

► Analysis: plotting P/Z = f{Gp) 

• The straight line becomes exponential in case of an active aquifer 

• The estimate of gas accumulation by extrapolation of the line is no longer valid => 
actually, at the beginning of the reservoir life, the influence of the aquifer may still be 
weak and the extrapolation may be valid yet uncertain => need monitoring (observation 
wells, 4D seismic) 




Gas reservoirs 

Havlenah-Odeh for gas reservoirs - Cole plot 

► The generalized MBE for a gas reservoir is given by: 
F = GE g + W e or F/E g = G + W e /E g 



where F = G p B g + W p B w and E g = B (J — B gi 

=> by plotting(F - W e )/E g = f(G p ) we may have a straight line of slope 0 and 
intercept G i.e. of equation y = G 


(F-W e )/E g 


G 




volumetric depletion 


strong aquifer 


moderate aquifer 


weak aquifer 


Cole plot is equivalent to Campbell plot for oil 
reservoirs. 


Note that E g should be replaced by E t if water and 
pore volume compressibility are not neglected: 

Et = Eg + Ef >w 

With 


Ef,w — Bgi 


^S wc ' c w 

i /r 


I A P 


G 
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Gas reservoirs 



Performance of a gas reservoir in case of water drive 

► In the case of natural water drive, the aim is to accelerate production in order 
to evacuate gas before water invades the reservoir and traps gas behind the 
front 

• Typical mobility ratio for water-gas displacement: 


M = 


krwmaw rgmaw 0- 2 


Pw 


Pg 


0 . 4 / 0 . 0.2 


= 0.01 


• This value means that, for a given A P in the reservoir, the gas travels a hundred times 
faster than water => favorable case but however some gas may remain trapped 
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Assessment of the trapped gas volume 
► Trapped gas volume 

^ trapped ~ BV inv 

With PV inv the pore volume invaded by water 


f grw 


l/B 


a 


And S grw the residual gas saturation to water that depends on the properties of the porous 


media/fluids system, especially wettability and IFT 

► By definition 


W — W B = PV 

e yy pMJ w i y inv 

= a ■ MGV 


(l $wc Sgrw^ a ■ P v(\ S wc Sg r w) 


With PV = GB gi /(l — S wc ) the pore volume 

MGV = PV(l - S wc - S grw ) the movable gas volume 
a the fractional volumetric sweep 
Therefore 


trapped 


W e - W p B w 

1 — s — s 

A J wc *“> grw 


f grw 


1 /B g = a- G 


' grw 


Bgi 


1 SwC Bg 
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Gas reservoirs 

Performance at abandonment pressure 



► Production at abandonment pressure is given 
by the intersect of the material balance line 
and the P/Z plot 


Some of the gas will be trapped at S rgw behind the 
front: 


aG Bgi S ° rw 

Bgab 1 


Some of the gas will be by-passed by the water: 


(1 - d)G 


} gi 


B 


gab 


Finally: 

f p,ab 


G„„ h = G- - (l - a)G^- 


Bgab 1 $wc 


B 


gab 


P 

%ab 



Zj f £grw_,lz£L\ 
Vl-S lvc + a ) 
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Gas reservoirs 

Gas condensates reservoirs 



► Dry/Wet gas reservoirs 




► Condensates gas reservoir 




Above dew point pressure, the previous equation is still valid 

Below dew point pressure, the liquids condensate in the reservoir and the previous 

equation has to be modified 



Pj % 2 phases 


Where Z 2 phases ' s th e two-phase compressibility factor taking into account the presence of 

liquid in the reservoir ; it is determined through CVD experiment 

And G p ' is the total production of hydrocarbons expressed in vapor phase: 




+ N 


p cond vap 
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Gas reservoirs 

Gas condensates reservoirs - 2 


► Total hydrocarbon production 

G„ — G n 4- N 


•p 


•p 


p cond vap 

23.96 

Np cond vap ~ ^ pPo * 


M 


w 


Where M w is the molar weight: M w = 
And °API = — 131. 5 

d 


6084 
°API- 5.9 
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Vf g = gas volume in R.C 

V gs = gas volume at S.C. 

V os = condensates volume at S.C 

V go - volume of condensates expressed in vapor 
phase (nb moles *mole volume std) at S.C 

Bg = Vfg/Vgs Bg = V fg/(Y gs + ^ go ) 
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Gas reservoirs 

Gas condensates reservoirs - Condensate banking 

> Condensates banking 

• Pressure drops mainly at the vicinity of the well 

• Hence, condensates will form first around the well where the pressure will first decrease 
below the dew point pressure => condensate banking / blocking 

> Possibly three zones in the reservoir 

• Away from the producing wells, monophasic gas 

• Close to the producing wells: condensate build-up region but with monophasic flow still. 
Liquids drop out of gas but remain immobile because their saturation is below the 
critical saturation S cc ; it is supposed to be a constant volume depletion region. 

• Closest to the producing wells, both phases are flowing since the liquids saturation is 
above the critical saturation 5 CC ; it is supposed to be a constant composition depletion 
region 

• This region ranges from a few tens of feet for lean gas to hundreds of feet for rich gas 

• The size is proportional to the volume of gas drained and to the percentage of liquid drop out 
(CGR) 

• It also depends on permeability => more permeability leads to more gas produced, and leads to 
more extent of the zone 

• Capillarity may play a role by retaining condensates up to a high values of saturation 
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Gas reservoirs 

Gas condensates reservoirs - Condensate banking 

► Managing condensates banking 

• Liquids should be of higher value than gas => the objective is to avoid condensates 
banking since part of the liquid will remain in the reservoir 

• Furthermore, condensate blocking disturb the production of gas 

► Solutions 

• Gas cycling: maintaining pressure over dew point by injecting gas, N 2 or C0 2 or dry gas 

• Finally dry gas (lean gas) may be produced at low bottom hole pressure 

• Economics: depends on the composition of the fluid and the quantity of liquids, on the 
availability of pipes to export the gas and on the liquid/gas prices 

• May be a seasonal activity or limited in time! 

• Limiting pressure draw-down in the reservoir 

• Increasing the contact area with the formation to increase production 

• Fracturing the reservoir (siliclastics) or acidizing the reservoir (carbonates) but condensates may 
form rapidly around fractures when sandface pressure drops below dew point pressure 

• Drill deviated/horizontal wells: longer time for the condensates to build up 

• Remobilizing the liquids 

• By injecting solvents but costly (product, monitoring, etc.) 
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Gas reservoirs 


0^ \ 

► Dry gas/wet gas reservoir 

• Without water entry and neglecting water and pore volume compressibility 

• Reservoir is produced through gas expansion only 

• Rf is a linear function of P/Z and can reach very high values, typically 90-95% 

• With water entry 

• Despite a favorable mobility ratio, gas may be trapped behind the water front 

• Rf is still a function of P/Z but no longer linear and will be lower than without water entry, 
typically down to 50-70% 

► Gas condensates reservoir 

• Above dew point pressure, it behaves like a dry gas/wet gas reservoir 

• Below dew point pressure 

• Need to use 2-phase compressibility factor to take the liquid into account in the reservoir 

• Rf will be lower for both liquid and gas because of the diphasic flow and the decrease in k r 

• Additionally, condensates banking: liquids will not flow until the critical saturation is 
reached 

• One solution may be gas recycling to maintain the pressure above dew point pressure 


IFPTraining 117 

Drive Mechanisms - Primary Recovery — ' 





Monophasic fluid expansion and pore volume reduction 




► It occurs when one phase only is mobile in the reservoir: 

• Undersaturated oil reservoirs as well as gas and gas-condensate reservoirs (above 
dew point) produce by fluid expansion and pore volume reduction. 


► Reservoir performance: 

• In undersaturated oil reservoirs, pressure declines very quickly; GOR remains 
constant until reservoir pressure falls below bubble-point pressure whenever GOR 
increases sharply (cf. solution gas drive); Rf is very low, typically a few % 

• In gas reservoirs, recoveries reach high value due to the combination of two specific 
gas characteristics, low viscosity and high compressibility; Rf is very high, typically 
ranging 90-95% 


V 
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Oil reservoirs - Solution gas drive 


\ 

► The solution gas drive mechanism describes oil displacement by the 

expansion of gas released from solution as pressure is reduced below 

bubble point: 

• As oil and gas production goes on, pressure declines further, more gas is released 
from solution. 

• Gas flow increases, oil flow decreases, as a result of the increasing gas saturation in 
the reservoir and of the unfavorable relative permeability evolution. 

► Reservoir performance: 

• Solution gas drive reservoirs exhibit typically strong reservoir pressure decline and 
correspondingly rapid oil production decline. 

• GOR rises quickly from initial value to a maximum value, before declining rapidly. 

• Oil recovery is generally rather low, typically ranging from 5 to 25 % 
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Oil reservoirs - Gas cap drive 


0^ 

► In presence of a gas cap above an oil zone, the pressure decline associated 

with the oil production will allow the gas cap to expand and provide 

additional energy to increase oil production 

• To be effective, a large gas cap is necessary (initial or secondary gas cap). 

• For a secondary gas cap to form, high vertical permeabilities associated with a 
relatively homogeneous reservoir are needed. 

► Reservoir performance: 

• An efficient gas cap drive mechanism exhibits typically a slower reservoir pressure 
decline and correspondingly a slower oil production decline than solution gas drive; 

• GOR rises slowly and progressively. 

• Oil recovery is generally higher than for solution gas drive reservoirs and it highly 
depends upon the vertical permeabilities. 

• Rf typically ranges from 20 to 40% 
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Oil reservoirs - Natural water drive 


CXu 

\ 

► For a water drive reservoir, the pressure primary source of energy is 
supplied by water influx (from an adjacent aquifer) into the reservoir. 

• In most cases, the energy for this water movement comes from water volume 
expansion and from the pore volume shrinkage in the aquifer. 

► Reservoir performance: 

• Water drive effectiveness is a function of the aquifer properties and not the reservoir 
ones. 

• The two key parameters are the aquifer size and the aquifer transmissibility 

• The total fluid rate may remain constant, and if reservoir pressure is kept above 
bubble-point, GOR remains constant. On the other hand, the reservoir shows a 
steady increase in WOR. 

• A model to estimate the water entries is needed, especially for large aquifers (e.g. 

Hurst and Van Everdingen model) 

• Oil recovery is high, typically ranging from 30 to 60 % 

V 
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► For a gas reservoir, the gas compressibility becomes dominant and is a very 
significant drive mechanism. 


► Reservoir performance: 

• In the absence of an active aquifer, high recoveries can be achieved, up to 90 %. 

• In presence of an active aquifer, the risk of trapping gas lead to lower recoveries, Rf 
may be as low as 60-70% 

• Produce as fast as possible 

• Observation wells needed to monitor the water rise. 

• Evaluate S grw (log - core) 


V I 
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Material Balance 



CXu 


fl 


► The material balance is a basic tool for reservoir engineers to: 

• Check production data consistency 

• Check the consistency between the geological evaluation and the reservoir behavior 

• Understand the reservoir behavior 

• Design a future complementary development plan 

• Perform production forecast 


V 
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Water influx models 


Schilthuis steady-state model 

► Schilthuis model assumes that, under steady state conditions, the aquifer 
flow can be described by Darcy's equation 

► Hence the rate of water influx can be expressed as: 


e w = dW e /dt = 


0.00708 kh 

H w In (r e /r 0 ) 


( Pi~p) = C ■ (Pi - p) 


where e w is the rate of water influx in bbl/d 
k is the aquifer permeability in mD 
h is the aquifer thickness in ft 
r e is the aquifer radius in ft 
r G is the reservoir radius in ft 

C is the water influx constant and can be calculated from the production history 
data provided that e w is calculated independently 
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Water influx models 



Schilthuis steady-state model - Solution 

► From previous equation, it can be seen that 

• W e = C Jfai - p)dt 

• The solution can be calculated by integrating p t — p as a function of time, for instance 
through graphical methods (area calculation) or any other integral approximation 
(trapezoid) 

• f‘(P,-p)dt = ^ t 1+ ■ (t2 - tl) + (Ecsslto). (tj _ tz) + 


► The cumulative water influx can be approximated by: 

W e = Ap • At 
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Water influx models 



Hurst modified steady-state model 


► One limitation of Schilthuis model comes from the aquifer radius 

• Indeed, the aquifer drainage radius increases with time 

• Hurst proposed to introduce a time-dependent dimensionless radius 


r e D = r e /r 0 = at 


• The corresponding modified water influx rate equation writes: 

0. 00708 kh ( Pi — P ) 

e w = dW e /dt = (Pi-p) = C 


Therefore 


p w ln(at) 

['(Pi 

Jo ln 


ln(at) 


W P = C I dt 


ln(at) 


Which can be approximated as: 




L 

"‘If 


Ap 


ln(at). 


At 
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Water influx models 

Solving Hurst modified steady-state model 

► Hurst modified steady-state model has two unknown values: C and a 

► They are estimated from reservoir/aquifer pressure history and from water 
influx rate history: 


Pi — pi,' 1 1 

= — In (at) = — In a + — In t 

e w C C C 


• Plotting ( pi — p)/e w as a function of In t should give a straight line of slope 1/ C and 
intercept 1/C ■ Ina 

• However, this method requires to get independent estimates of e w as a function of time 
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Water influx models 



HVE model - Solution for edge-water drive linear aquifer 


• Aquifer constant 
U = wLh0c t Darcy units (cc/atm) 



U = 0. 1781wLh0c t Field units (bbl/psi) 



Dimensionless time 


Darcy units (t-seconds) 


Field units 


with c te = 0. 000264 t-hours 
c te = 0. 00634 t-days 
c te = 2. 309 t-years 


Note that t D is now a function of the size of the aquifer 
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Water influx models 

HVE model - Bounded aquifers 



► For bounded aquifers, it can be shown that W D reaches a constant maximum 
value 


• Radial aquifer: W eD (max) = 0. 5 {r eD 2 — l) 

In this case and for a 360° aquifer, it can be shown that: 

W e = n (r e 2 — r o 2 )h0c t Ap = W t c t Ap 

We find back the small pot model i.e. the total influx of water assuming the pressure drop 
A p has been transmitted instantaneously throughout the aquifer 

• Linear aquifer: W eD (max ) = 1 
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Water influx models 

HVE model - Infinite aquifers 



► No maximum value of W eD can be reached since the water influx is always 
under transient conditions 

• Radial aquifers with water-edge drive: Edwardson & al. (1962) developed 3 sets of 
simple polynomial expressions for calculating W eD : 

• Fort D <0.01 W c d = V^d/ 71- 


• For 0.01 < t D < 200 


W eD = 


1.2838Vf^+1.9328t D +O.269872t0 3/2 +O.OO855294t g 2 

l+0.616599Vt^+0.0413008t D 


• For t D > 200 W eD = (-4. 29881 + 2. 02566t D )/ln t D 


• Linear aquifer: it is possible to calculate the water influx directly 
W e = 2hWy[0kc^t/np x A p in Darcy units (cc/s) 

W e = O.OO326huV0^7^ x ' n un ' ts (bbl) 
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Water influx models 

Fetkovitch method - Approximated HVE model 
► Principles: 

• Introduced in 1971 in order to simplify the calculation of HVE model 

• The flow of water from the aquifer into the reservoir is modeled the same way as the 
flow of oil from the reservoir into the well: 

dW e 

q w = ~^r = KVa ~ P) 

where: q w is the water influx rate 

/ is the aquifer productivity index 

p is the reservoir pressure i.e. the pressure at the WOC (or GWC) 
p a is the aquifer average pressure 

• Average pressure in the aquifer is defined as: 

W e cV\v(Pi Pa) 

where p f is the initial pressure in the aquifer AND the reservoir that is 

with W ei = c t WiPi is the initial amount of encroachable water 

/PP Training I 134 

Drive Mechonisms - Primary Recovery — ' 


Water drive 

Fetkovitch method - Solution 

► After various steps of differentiation/integration we get: 

W e = — — (pi — p) ■ (l — e~ JPlt ^ Wel ) 

Pi 

When t -» 00 then W e -> — - (p t - p) = c t Wi(pi — p) 

Pi 

That is the maximum value of water influx that can occur once the pressure drop p f — p has 
been transmitted to the whole aquifer 

> Practical use of Fetkovitch method 

• The previous solution is only valid for a constant inner boundary pressure and we still 
shall need to use superposition theorem in order to use it 

• However Fetkovitch showed that superposition could be avoided using a different form 
of the previous equation 
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Water influx models 

Fetkovitch method - Practical solving method 
► Methodology 

• During the first time step A t ± we get: 

(Pi - Pi) ■ (1 - 

where p 1 is the average reservoir boundary pressure during the first time step 

• During the second time step A t 2 we get: 

A W e2 = SjS (p ai - p 2 ) ■ (1 - e~IP^/w el) 

where p aj = p f ( 1 ) is the average aquifer pressure during the first time step 

^ ™ ei ' 

and p 2 is the average reservoir boundary pressure during the 2 nd time step 

• In general, for the n th time period we get: 

AW en = ^ (p an _ t ~ p n ) • (l - 

1 

where p anl = P/ 1 1 — — - I is the average aquifer pressure during the n-lth time 

step 

and p n is the average reservoir boundary pressure during the nth time step 
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Water influx models 

Fetkovitch method - Aquifer productivity index 



► In order to be solved, Fetkovitch method requires that the aquifer 
productivity index / is calculated 

► Pseudo-steady state with A P = Vi~Pa 

• Radial aquifer 

J = 7 r ~ — \ Darcy units (cc/s/atm) 

K ln ?r 0 - 75 ) 

J = 0. 00708 , { kh — t Field units (b/d/psi) 

*( ln f 0 -°- 75 ) 

• Linear aquifer 

/ = 3 Darcy units (cc/s/atm) 

/iL 

J = 3.381^y- Field units (b/d/psi) 

flL 
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Water influx models 
Fetkovitch method - Limitations 


► Small aquifers 

• The expression of aquifer productivity index / used in Fetkovitch modeling is valid only if 
(r 0 /r e ) 2 is negligible 

=> this can be untrue for small radial aquifers but Fetkovich approximate showed robust 
behavior to this assumption 

► Large aquifers 

• The initial influx of water is under transient flow conditions and the corresponding 
expression of aquifer productivity index / is no longer valid 

=> Fetkovich method cannot model infinite aquifers 

=> for very large aquifers, it is still necessary to use the FIVE method at early times 
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Drive mechanisms 
Secondary recovery 
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Secondary 

Recovery 


Introduction to secondary recovery 
Principles 



► The natural drainage of oil reservoirs often leads to low recovery factors 

• Additional energy has to be injected into oil reservoirs in order to achieve higher 
recovery 

• First solution is to inject fluid, water and/or gas, in the reservoir => Secondary recovery 


► Objectives: to increase the reserves 

• By supporting the pressure 

• By improving the sweeping efficiency of the hydrocarbon zone 

► Methods 

• Water injection (in the aquifer or in the lower part of the oil zone) => water flooding 

• Gas injection (in the gas cap or in the top on the oil zone) => gas flooding 
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Introduction to secondary recovery 
Example of secondary recovery 




From Sneider 


Drive Mechanisms - Secondary Recovery 


!PP Training 


14 

3 


Introduction to drive mechanisms 



Steam injection 
In situ combustion 
Electric heating 
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Introduction to secondary recovery 
Water injection / Gas injection 



Producers 


► Water injection: 



Introduction to secondary recovery 
Evaluation of drive mechanisms 



► Reservoir engineers have to determine from the onset the reservoir production 
mechanisms: natural or assisted drainage (water/gas/steam injection,...) 

► Generally, assisted drainage is preferred because of higher recovery factors but 
natural drainage mechanisms always have to be determined in order to 
optimize assisted drainage 

► It is interesting to start producing the reservoir through natural depletion, even 
for a very short period, in order to observe its behavior and determine the type 
of natural drainage mechanisms underway from the analysis of production data 

► However, assisted drainage is more and more often implemented right from 
the beginning in order to achieve a higher recovery, 

• Especially since early recovery is highly valuable 

► As engineers learn more about the reservoir, the initial FDP may be changed 

• Adding new wells, either producers or injectors, converting producers to injectors, etc, 

► Nevertheless, uncertainties may remain, in particular aquifer activity 
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Introduction to secondary recovery 
Parameters to define /optimize 

► Many parameters have to be defined / optimized: 

• Nature of the fluid to be injected (water, gas, steam, ...) 

• Injection zone (top, bottom, aquifer, all...) 

• Injection pattern (grouped flood, dispersed flood,...) 

• Optimum level of pressure support 

• injection rate, continuous/alternate injection,... 

• number of injection wells 

• planning 

• Field surveillance and data acquisition to monitor the results of injection strategy 
implementation vs. forecast 

• Data input for well architecture/completion design and for surface facilities design 
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Introduction to secondary recovery 




► The main objective of secondary recovery is to increase hydrocarbon 
reserves by: 

• Supporting pressure in the reservoir 

• Improving the sweep efficiency of the reservoir 

► The most commonly used solution is to inject fluid: 

• Water injection, in the aquifer or in the oil zone near the WOC 

• Gas injection, in the gas cap or in the oil zone near the GOC 

► Secondary recovery implies multiphase flow in the reservoir 

• Water/oil displacement in case of water injection 

• Gas/oil displacement in case of gas injection 

► To be optimized, assisted drainage requires information on natural 
drainage, but it is often implemented without sufficient appropriate 
information 

• FDP will be updated all along the reservoir life to optimize the secondary recovery by 

V taking information coming from dynamic/production data analysis into account 
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Oil reservoirs drive mechanisms and injecting fluids 


► Which kind of drive mechanisms are favorable for secondary recovery? 

• Undersaturated depletion 

• Waterflooding: it may prevent pressure to drop below bubble point and avoid free gas flowing 
in the reservoir; it will also improve Rf otherwise very low 

• In some cases, allowing pressure to drop below the bubble point can be favorable because of 
triphasic effects (solution gas mobilizing the oil and slowing down the water) 

• Solution gas drive reservoirs 

• Very good candidates for waterflooding since these reservoirs may exhibit rapid pressure and 
production decline and still quite low recovery factors 

• May also be good candidates for gas injection because of the large quantity of associated gas 
available at the beginning of the process 

• Gas cap reservoirs 

• Good candidates for waterflooding but attention has to be paid to avoid losing water in the gas 
cap or pushing the oil into the gas cap 

• May be also very good candidates for gas injection since associated gas is available and it is 
possible to inject directly in the gas cap with very good injectivity 

• Natural water drive reservoirs 

• Waterflooding may be applicable in the case of weak water drive in order to improve pressure 
support and recovery 
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Secondary recovery projects 
Nature of injected fluids 

► Performance of injection 

• One of the main driving parameters of the fluid injection performance is the mobility 
ratio M that has to be as low as possible in order to avoid instabilities hence a quite 
viscous injection fluid 

• As a consequence water is well adapted for light oils, but less for heavy oils 

• In any case, injected gas poorly sweeps oil in place due its lower viscosity but immiscible 
gas injection can be turned into miscible gas injection and achieve better performance 

• Generally water injection is performed first, then gas injection, before moving possibly 
to EOR in order to further increase the performance (miscible gas, WAG and finally 
possibly chemical e.g. SP, ASP, FAWAG depending on the case etc.) 

► Treatment 

• Water treatment is generally easier than gas treatment 

► Availability 

• Key factor for implementing a secondary recovery injection scheme 

• Water availability is generally higher than gas availability (mainly non-marketed 
associated gas) but gas injection can be seen as a way to increase value after the 
prohibition of venting or flaring 
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Secondary recovery projects 
Injection pressure 

► One of the main objectives of water injection is to maintain (or support) the 
reservoir pressure, especially above bubble point pressure in order to avoid 
free gas in the reservoir 

• To prevent the increase in oil viscosity by « losing » the lighter components 

• To prevent the decrease in oil relative permeability because of the increase in gas 
saturation 

• To maintain good wells productivity (IPR) 

• To prevent a sharp increase in GOR and related problems at the surface 

► As a rule of thumbs, the typical target for pressure maintenance should be a 
few bars to a few tens of bars above bubble point pressure 

• High pressure maintenance requires high pressure head and may be costly 

• However, it is much easier to inject water than gas because of the higher value of B w 
compared to B g that decreases roughly as 1/P => as a consequence, the volumes of 
water to deal with at the surface are much lower and the corresponding equipment is 
less costly 

► The target in pressure may typically rule the start of the injection period 
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Secondary recovery projects 
Volume of fluid injected 



► For the reservoir pressure to be strictly maintained, there should be a strict 
balance between the injected fluid volume and the fluid voidage at reservoir 
conditions 

► If no other mechanism is involved (e.g. active aquifer), the injected fluid 
volume should compensate the produced volumes at reservoir conditions: 


injection rate = production rate @ res. conditions 
Qwinj^winj + QginjB ginj ~ Qo^o “h Qwp^wp ”1” QfgBfg 
where Qf g is the surface flowrate of free gas Qf g = Q 0 (fiOR — R s ) 


► Actually, the injected volume will also depend on the availability of the fluid 
for injection 
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Secondary recovery projects 
Voidage replacement ratio 



► It is the ratio of the injected volumes to produced volumes in reservoir 
conditions 

• Instantaneous VRR 

QwinjBwinj “h QginjB ginj 


VRR = 


QoBo + Q o (G 0R — R s )B g + Q W B W 


Cumulative VRR 


VRR = 


W injB winj "h ^injBginj 

N p B 0 + N p (R p — Rs)B g + WpB w 


• VRR target should be 1 or above 1 in order to maintain pressure 

• In some cases it can be interesting to let the VRR fall well below 1 in order to let the pressure 
drop below bubble point and benefit from the triphasic effects 

• Cumulative VRR are typically calculated on a monthly or three-month basis and since the 

beginning of the injection process 

• Monthly cumulative VRR may help to identify injection losses or unexpected fluids influx into 
the reservoir 

• As long as the overall cumulative VRR is equal to or above 1, it can be said that the reservoir 

pressure has been maintained from the beginning of the injection process iFPTraining 
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Secondary recovery projects 

Voidage replacement ratio and oil flowrate 



► Plotting VRR together with the oil rate may help understand the relationship 
between the two 

• The below figure displays VRR and the oil flowrate (log scale) vs. time for El Trapial field 
(Argentina) and shows a clear relationship between the two 
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Secondary recovery projects 
Injection pattern 



► The injection pattern, i.e. the relative layout of injection and production wells 
is critical for the success of a waterflooding / gas flooding project 


► It strongly depends on the reservoir geometry, the reservoir heterogeneity 
especially in permeability, the fluid content and the volume to be swept 


► There are two kinds of injection patterns 

• Grouped flood where injection wells are grouped locally 

• Dispersed flood where injection and production wells are alternatively located 
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Secondary recovery projects 
Grouped injection pattern 



► In the case of high-dip reservoir (gravity effects) or high permeability 

reservoir (with associated low pressure gradient), injection wells are located 
to get a regular displacement i.e. piston-like allowing to sweep large areas 
with late breakthrough at the production wells => grouped flood 


• In case of gas, injection is made 
close to the GOC or in the gas-cap 
if any, or at the top of the 
structure => central flood 



• In case of water, injection is made 
in the aquifer or near the WOC => 
peripheral flood 
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Secondary recovery projects 
Dispersed injection pattern 



► If the reservoir shows no gravity effect (horizontal) or exhibits low 
permeability, only a very limited zone is effectively flooded 

=> production and injection wells are laid out in a fairly regular pattern in the oil zone 
termed dispersed flood 


► Several patterns are considered as characterized by their geometry and the 
ratio injectors/producers: 

• I / P ratio 

• 5 - spot I = P 

• 7 -spot I / P = 2 regular 

• I / P = 1 / 2 inverse 

• 9 - spot I / P = 1 / 3 
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Secondary recovery projects 

Dispersed injection pattern - Wells patterns 
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Secondary recovery projects 
Wells locations 



► 


► 


► 


The production/injection scheme 
is designed according to the 
reservoir characteristics 

• Distance between wells 

• Relative location of injectors and 
producers 

Due to high complexity 
(heterogeneities, faults, etc.), 
reservoir simulation is generally 
needed to determine actual 
optimal well locations 

With time: 

• Producers may be turned into 
injectors 

• Additional injectors may be drilled: 

infill drilling 



Drive Mechanisms - Secondary Recovery 


!FP Training 


160 


Secondary recovery projects 
Waterflooding - Some history 

► Waterflooding is the most widely used secondary recovery technique 

► Water injection began in the U.S. as early as the end of 19th century and it 

developed widely early in the 20th century 

• The first concern was about production water management: as the fields became 
mature, more and more water was produced and disposal in rivers became more and 
more difficult => start of water injection included within the reservoir interval 

• Injection of water was generalized by the 1920-1930s 

• In parallel, engineers recognized the very low recovery factors obtained from primary 
recovery and the improvement of recovery due to water injection 

• By the 1940s, secondary recovery with water injection began to be systematically 
implemented in the U.S. and engineers began to work on modeling and optimizing water 
injection 
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Secondary recovery projects 
Waterflooding - Injectivity 



► Well injectivity is one of the major concern and source of uncertainty within a 
secondary recovery project 

► Well injectivity depends on the location: 

• In the water zone, Injectivity Index II is given by: 


II = 


W 


P . . — P 
r win] r r 


= 0.00708* 


kh 




w 


ln^ + S 

r w 


In the oil pool, Injectivity Index II is given by: 


II = 


:w 


P — P 

1 win j 1 r 


= 0.00708* 


IdCyyy h 




w 


ln^ + S 

r w 


=> The injectivity index is better in the water zone since k rw « 0. 3 but if injecting in the 
water zone, the effectiveness will strongly depend on the quality of communication 
between the aquifer and the oil pool (especially vertical permeability) and the quality of the 
permeability distribution in the water zone (largely unknown) 
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Secondary recovery projects 
Waterflooding - Practical considerations 

► Water sources 

• Fresh surface water 

• Offshore: sea water 

• Water from a different reservoir (above, below, ....) 

• Produced water 

► Potential problems 

• Compatibility with formation water (BaS0 4 , CaC0 3 precipitates) 

• Water filtration 

• Biocides to prevent H 2 S generation 

• Oxygen removal 

► Injection of tracers for monitoring the waterflooding 
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Secondary recovery projects 

Waterflooding - Performance and related parameters 

► Volumetric recovery 

• Recovery at the reservoir scale, by respect of the part of the volume that is contacted 

• Derived from the flow simulation (cross section & full field model) 

• Can be high in the case of a favorable mobility ratio M 

► Microscopic recovery 

• Recovery at the scale of the pore, i.e. the part of the pore volume that is filled with 
water within the pores which are contacted by the waterflood 

• Derived from the relative permeability curves 

• Generally quite low because of high value of residual oil saturation 

► Overall performance 

• Can be as high as 50-60% 

• Driving parameters 

• Permeability distribution 

• Relative permeability & water/oil mobility ratio 

• Faults - Compartmentalizing 

• Replacement rate 
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Secondary recovery projects 
Waterflooding - Gravity drainage 

► Other mechanisms may have an impact on recovery, depending on the 
context (geology, fluid): 

• Capillarity 

• Gravity drainage 

► Gravity drainage can dramatically enhance the sweep efficiency by stabilizing 
the flow 

► Gravity drainage is favored by the following situations: 

• High reservoir dippage 

• Strong vertical permeability k v and no vertical permeability barriers 

• High water-oil density contrast 
• Better with light oil 
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Secondary recovery projects 
Gas injection - Performance 

► Main drawback 

• Due to low viscosity ( n g = 0. 01 cP ), the mobility ratio M is high, instabilities develop 
and the volumetric sweep efficiency is low => only a small part of the area under 
injection is effectively flooded 

• Highly sensitive to permeability heterogeneities 

► The overall performance is generally lower than water injection but 

• A high A p og leads to good segregation especially in case of reservoir dippage and, if k v 
is high, gas injection results in a very efficient gravity drainage with a recovery as high as 
60-70% 

• A higher gas mobility gives capacity to flood a part of the reservoir that has not yet been 
swept by the water 

• The residual oil saturation to gas S org can reach very low value, in the range of 5%, thus 
leading to high microscopic efficiency 

• Additionally immiscible gas injection can be turned to miscible gas injection achieving very high 
microscopic efficiency 

• Gas injection may come in support of waterflooding in order to further increase the 
recovery by remobilizing isolated oil and flooding areas kept unswept 
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Secondary recovery projects 
Gas injection - Practical aspects 

► Injectivity 

• Taking into account low ii g , injection wells have high injectivity => less gas injectors 
(typically less than for producers) 

• But necessity of high wellhead injection pressure => compressors 

► Availability 

• One of the main concerns: finding gas 

• Large volumes of gas needed 

• Most of the time, the injection gas comes from produced gas (associated gas) 

• Associated gas was usually flared but flaring has become more and more difficult 
because of environmental concerns => reinjection can be a valuable solution but the 
increasing concern for "marketing" associated gas may become a problem for secondary 
recovery project 

• Injection may also be a solution to valorize gas in the case of gas condensates reservoir 
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Secondary recovery projects 


0^ 

► Waterflooding is the most widely used secondary recovery method since 

• The availability is generally better 

• The treatment is generally easier 

• The sweeping efficiency is generally better 

• The overall performance can be as high as 50-60% 

► Gas injection may be used to further improve recovery 

• The sweeping efficiency is generally lower but microscopic efficiency is better 
because of lower residual oil saturation to gas 

• A higher mobility allows to flood areas unswept by water and may help to remobilize 
isolated oil 

• The overall performance can be as high as 60-70% in case of a favorable gravity 
drainage, otherwise generally lower than for waterflooding 

• Immiscible gas injection can be turned to miscible gas injection with even better 
performance (cf. EOR) 

• Main concern: gas availability and related cost 

V 
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Secondary recovery projects 

CXu 

\ 

► The injection pattern is of two types: 

• Grouped flood in the case of high-dip reservoir or high permeability reservoir: central 
flood in case of gas injection, peripheral flood in case of water injection 

• Dispersed flood in the other case using patterns like 5-spot, 7-spot etc. 

• In general, a reservoir simulation is used to optimize well locations 

► Well injectivity 

• Critical but very often difficult to forecast 

► Pressure 

• Typically a few bars to tens of bars above bubble point pressure 

► Overall performance strongly depending on the reservoir geology 

• Deep knowledge of the reservoir mandatory, especially heterogeneity 

• Monitoring of the performance over time (cf. VRR) 

► Surface facilities, especially treatment facilities, have to be designed 
accordingly 
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Mobility ratio 

Mobility ratio and flow stability 

► Instabilities in multiphase flows 

• When the reservoir is thick and vertical velocities cannot be ignored, it can be shown 
through the analysis of forces that distortions or encroachments may appear and 
develop at the fluid interfaces 

• These encroachments occur at different scales: 

• fingering at small scale 

• tongues at the scale of the reservoir 

• coning near the well 

• Encroachments are governed by conditions of stability or instability of the flow 

• Among the parameters governing instabilities is the mobility ratio M: the lower the 
mobility ratio, the better the stability of the displacement 

• The lower the mobility contrast between the displacing fluid and the displaced fluid, the more 
stable and the more efficient the displacement 
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Mobility ratio 

Favorable and unfavorable mobility ratio 



► It can be shown in laboratory that fingering starts to develop at microscopic 
scale when M > 1 and tends to disappear when M < 1 


► A mobility ratio M above 1 is said to be unfavorable while a mobility ratio 
below 1 is said to be favorable 
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Mobility ratio 

Calculation of mobility ratio 



► Definition 


In the case of oil/water displacement: 

k rw C^vvra ) Mo 

M = . x 

Mw K'ro Wwi) 


where is the average water saturation behind the front 
and k rw (S wm ) is the corresponding water relative 
permeability (from k r curve) 

• If the viscosity contrast between the two fluids is not too 
high, then S wm is close to S wmax = 1 - S orw (i.e. the 
endpoint of the k r curve) => end-points mobility ratio 

• Otherwise, we should need to use S wm estimated from 
Welge tangent analysis (see below) 

• Actually, the end-point mobility ratio is generally used to 
decide on the stability condition 


01 
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Relative permeabilities 

K r curves - Example of oil/water system 




1: S w decreases => drainage (corresponding to the 
initial flooding of oil in the reservoir supposed 
initially water-wet) 

2: S wirr irreducible water saturation : only oil moves 
k ro =0.8 & k rw =0 

3: S w increases => imbibition (corresponding to oil 
production) 

4: S or residual oil saturation : only water moves, oil 
is immobile 

5: crossing point : k ro =k rw and k ro +k rw < 1 

6: hysteresis : not the same curve for imbibition and 
drainage due to difference in invasion process 
between wetting fluid and non-wetting fluid 
(effect of capillary pressure); hysteresis is 
generally more marked for the non-wetting fluid 


• Except for the end points, the sum of the relative permeabilities is always < 1 => each 
fluid disturbs the flow of the other one 

• K r curves dictate the fluids flow in the reservoir and are used for reservoir simulation 
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Mobility ratio 

Examples of typical values of M 




From Cosse 
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Mobility Ratio 



<Xu 

'N 

► Mobility ratio M is defined as the ratio of the mobility of the displacing 
phase to the mobility of the displaced phase 

► M is one of the main parameters governing the instabilities which may 
develop during multiphase flow and affect the sweep efficiency (tongue, 
fingering) 

► If M < 1 it is favorable ; if M > 1 it is unfavorable 

• M is favorable only in some cases for light oil displaced by water; displacement of oil 
by gas always leads to an unfavorable gas ratio due to very low gas viscosity 

• An unfavorable M leads to a low value of S wm far from maximum water saturation 
=> more oil trapped behind the displacement front 

V 
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Water drive under 
diffusive flow 
conditions - 
Buckley-Leverett 
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Diffusive flow - Buckley-Leverett theory 
Water injection / Gas Injection 




► Gas injection 

• Into or close to the gas cap if 
any, at the top of the 
structure otherwise in order 
to form a secondary gas cap 

• The gas cap is getting down 
into the reservoir => one 
mobile interface gas-oil 
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► Water injection 

• into or close to the WOC 

• The WOC is rising into the reservoir 

=> one mobile interface water/oil 


Gas injectors 



Diffusive flow - Buckley-Leverett theory 
Main assumptions of diffusive flow 



► Immiscible fluids 

• Typically water displaces oil in a water-wet reservoir 


► Displacement occurs under condition of vertical equilibrium 

• Saturations are distributed under capillary-gravity equilibrium according to the 
equation: 

P C (S W ) = A p. g. cos 0. h 

• The vertical velocities are considered as infinite compared to the horizontal velocity 
=> S w is instantaneously redistributed by respect to the previous equation 

=> fluid saturations are uniformly distributed with respect to thickness 

=> unidirectional displacement of a water saturation front 


► Displacement is considered to be incompressible 

• Total fluid replacement => reservoir pressure is maintained 

fit Qo Qw Qitij 

► Displacement is linear 


Drive Mechanisms - Secondary Recovery 


!FP Training 


182 


Diffusive flow - Buckley-Leverett theory 
Conditions for vertical equilibrium 

► Vertical equilibrium is enhanced by 

• A high value of permeability, especially vertical permeability 

• A high value of density contrast between fluids 

• Low fluid viscosities 

• A low reservoir thickness 

• High value of capillary forces, i.e. large transition zone 

► Diffusive flow applies typically 

• In plug/core flows, in laboratory 

• Within the reservoir with large transition zone where there are major capillary forces 

► But was found valid with good results within most of the reservoirs 
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Diffusive flow - Buckley-Leverett theory 
Remainder on relative permeabilities 



► The classic Darcy's law is valid for monophasic flow only and it needs to be 
extended to take the reservoir multiphasic flows into account 


► Case of diphasic flow 

• Two fluids flowing in the porous medium 

• The monophasic Darcy's law is extended by introducing the effective permeability that 
takes the presence of the other fluid in the porous medium into account 


Qi = j a ^T where is the effective permeability of fluid 1 by respect to fluid 2 
Q 2 = where k 2 is the effective permeability of fluid 2 by respect to fluid 1 
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Diffusive flow - Buckley-Leverett theory 
Kr-Pc reminders 




Qw = “ 21 with k rw = k w /k 



0.0 0.2 0.4 0.6 0.8 1.0 


► Water saturation, S w ◄ 

s s 

wi ^or 
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Diffusive Flow - Buckley-Leverett theory 
Frontal Displacement 



► Case of two immiscible fluids (oil and water); unidirectional flow 



Drive Mechanisms - Secondary Recovery 


Zone I: Unreached by the displacing fluid: oil is the only fluid in movement 

Zone II: Front : S w increases rapidly 

Zone III: S w varies gradually behind the front 

Zone IV: Flooded by water: water is the only fluid in movement; residual oil S or 
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Diffusive flow - Buckley-Leverett theory 
Ideal and real displacement 


Ideal displacement: piston-like displacement 




1 



LJ 


DDE 

ISor 



III 

“5 

S W 

0 



0 x L 



Real displacement: function of k r , ju, k, p, q, dip 



Leaky Piston type 
displacement 
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Diffusive flow - Buckley-Leverett theory 
Unidirectional frontal displacement 


► One can demonstrate that within a section of 
width dx, the ratio q w / q t depends on the 
water saturation S w only 


► 


Darcy's law writes (without gravity): 


q 0 = - 
q w = - 


k 0 A dP 0 

Ho dx 
k w A dP w 

Pw dx 


► 


Gravity introduces an additional term along 
with the pressure gradient in Darcy's law 


k 0 A (dP 0 . > 

q ° = -j;{iu +p ° asina j 

k w A ( dP w 

qw = ~ + P'»9 slna 
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Diffusive flow - Buckley-Leverett theory 


Fractional flow equation 

► Introducing capillary pressure P c = P 0 — P w we get: 

-^ = -5f + P»g sma W 



and 

Qw\*w d(P q-P c ) 


+ p w gsina (2) 


k w A dx 

► Subtracting (1) from (2) we get: 

~dP c 


k w A k 0 A 


dx 


+ Apgsina 


q w + q 0 = - Apgsina ] (3) 


► Introducing fractional flow (or water-cut) f w in equation (3) 


fw = 


Qw 


— „ Qw fwQt ^nd (j a — Qf (1 f w) 


fw = 


Qw+Qo Qt 

. , Ak 0 (dP c A . \ 

1 + — —•(- r g —A pg sin a 
_ PoQt \ dx ™ ) 


1 + 


Pw/kp 

ppk w 


As P c , k Q and k w depend only on f w is a function of S w only 
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Diffusive flow - Buckley-Leverett theory 
Simplification of fractional flow equation 

► Simplification assuming only viscous forces: 

• Neglecting capillarity => P c = 0 

* Ak 0 A 

1 — A p.g .sin a + r 

f _ PpQt _ 1_Cj 

' w ^ | Pwkp i Pwkp 

Ppk w Ppk w 

with G the gravity term 

G = — - ■ A p. g . sin a 

PoQt 

• Neglecting gravity (horizontal reservoir or low dippage) 


' w ^ | Pwkrow(Sw) 

l*okrw(.Sw) 

• In this latter case, f w (S w ) is calculated directly from fc r (5 w ) tables => a very easy way 
to evaluate the fractional flow/ w as a function of S w 
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Diffusive flow - Buckley-Leverett theory 
Fractional flow curve 



► From the curves k r (S w ) — P C (S W ), we derive the curve f w (S w ) 
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Diffusive flow - Buckley-Leverett theory 
Buckley-Leverett equation 



► Considering a displacement of oil by water, the mass balance writes: 


q w 


■> 



Ax 


I 


) 


[( ^lwPw)x (.QwPw)x+dx\^t A0Ax\(S^ vPw) (S\vPw) ] 


^ ^ 

When Ax -> 0 and At -> 0 , we get: — ( q w p w ) = A0 — (S w p w ) 
Assuming p w — c ste and introducing q w — f w q t we get: — 

Since f w = / W (5 W ) we finally get the Buckley-Leverett equation: 

-df w dS w _ A0 dS w 
dS w dx q t dt 
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Diffusive flow - Buckley-Leverett theory 
Constant saturation frontal advance equation 



► Since S w = S w (x, t) and as we follow a fluid front of constant saturation, we 
get: 

dS w dS w 

dS w = —^dx + —^dt = 0 
w dx dt 

► Introducing this equation in the previous Buckley-Leverett equation, we get: 

dx _ q t df w _ . 

dt A0dS w 1 w) 

=> the velocity of a front of a given saturation is constant 


► Integrating in time, we get: 

_ fit ' t (df w \ _ Rinj ' t fdf w \ 
Xs -'~~m\ dS^,)' A® \dS^J 

where W in j is the volume of water injected so far 


Wjnj ( dfw\ 

A0 \dS w ) 
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Diffusive flow - Buckley-Leverett theory 
Buckley-Leverett solution 



► From f w (S w ) and its derivative, it is possible to get S w vs. distance x 

• The plot clearly displays an impossible physical situation with two different values of S w 
for a given x (due to the presence of an inflexion point in the fractional flow curve) 
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Computed water saturation profile 
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Diffusive flow - Buckley-Leverett theory 
Buckley-Leverett solution - 2 



► Buckley-Leverett solution consists in introducing a frontal discontinuity in £„,(;*:) 

• The areas behind and before the front are balanced 

• Physical interpretation is that all the water injected remains behind the front 


Balancing of areas 
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Diffusive flow - Buckley-Leverett theory 
Situation before and at water breakthrough 

Before breakthrough 



At breakthrough 
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Diffusive flow - Buckley-Leverett theory 
Water saturation at the front 



► In general, we can write 


r$wf 

%f(S w f — S w i ) — I x. dS w 

$wi 


► From the previous equation we have: 



► Furthermore 

_ Qt_ 't { df w \ 
40 \dS w J 


=> 


!sj X dS ™ ~ 40 fsj ~ 40 ^ w( S wf) fw( S wi )) 


► Therefore: 

(df w \ - /»,( 5 wi ) 

\dS w ) Sw/ $wf ~ Swi 

=> the slope of the tangent of the curve / W (5 W ) at the point S W f is equal to the slope of the 
chord joining the first point of the curve / W (5 W ) (where it is assumed that S w = S wi ) and 
this point => graphical construction using Welge tangent 
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Diffusive flow - Buckley-Leverett theory 

Water saturation at the front - Welge tangent method 



► Welge tangent is the tangent at the 
curve f w (S w ) passing by the initial 
water saturation point 

► The line passing by the point 
(.Sw = S wi ,f w (S wi ) = 0) and 

tangent to the curve defines the 
water saturation at the front 

► The Welge tangent allows to 
determine graphically the water 
saturation at the front (and 
corresponding fractional flow) 


Determination of saturation at the front 
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Diffusive flow - Buckley-Leverett theory 
Average water saturation after the front 



► The calculation of the slope of the 
Welge tangent gives: 

fdf w \ 1 

\dS w J S\vm $wi 

J wf 

where 5 wm is the value of the average 
water saturation behind the front 

► The intersection of the Welge 
tangent with the line f w = 1 gives 
the value of the average water 
saturation behind the front S wm 


Determination of the average saturation 
behind the front 
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Diffusive flow - Buckley-Leverett theory 


Oil recovery and time of water breakthrough 



► The oil recovery (in % of the pore volume) at breakthrough is given by: 

d/V 


Npd,bt Widbt S\vm,bt S\vi $ wm $ wi ^ 


dS 


Wj 


>wf 


and 


Rf = 


s — s 

J wm J wi 


1-5 


Wl 


► The velocity of the front is given by: 

y _ Qt Af w _ / df w \ (fa 

f ~A0AS^~A0\dS^) <: A0(S wtn - S wt ) 

3 wf 

► Given the distance L between the injector and producer wells , the time of 
breakthrough is: 

j _ L _ LA0(S win — Swi) _ LA0(S wm — S w i) _ W i n j _ Wi d 

Vf Qt Qinj flinj Olid 

where W id is the dimensionless injected water volume: W id = S wm — S wi 
and q id is the dimensionless injected water flowrate: q id = q in j/LA0 
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Diffusive flow - Buckley-Leverett theory 
Situation after water breakthrough 



► 

► 


At x = L, water saturation S we increases over S w f and fractional flow f we 
increases 

We caq still use Buckley-Leverett equation at x = L and we get: 


Wi (df w \ 
A0 \dS w ) 



0 L 
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Diffusive flow - Buckley-Leverett theory 
Oil recovery after water breakthrough 



► Given water saturation 
S we > S W f and fractional 
flow f we we may draw a 
Welge tangent and get the 
corresponding value of 
average water saturation 
S 

** wm 


► The oil recovery (in % of the pore volume) is given by: 

Npd S wm ~ S\vi (^wm ^we) (^we — $ wi ) (^we S wi ) d" (1 — fwe ) ' ^ id 

where W id is the % of pore volume invaded given by: W id = yy- = l/i-jf-) 

LA\D \ab w y r* 

3 we 




► The corresponding time (assuming constant injection rate) is: 

Qid Qinj ' \dS w J„ 
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Diffusive flow - Buckley-Leverett theory 
Effect of the mobility ratio 



► Lower favorable values of the 
mobility ratio M shift the 
fractional flow curve to the 
right 

• Swm value is high, almost equal 
to the maximum water 
saturation 

• Most of the oil has been 
recovered at the breakthrough 

► Higher unfavorable values of 
the mobility ratio M shift the 
curve to the left 

• Swm value is low, far from the 
maximum water saturation 

• Large quantities of oil remain 
behind the front 



Sw 


Reminder: simplified fractional flow equation is: 


f w 


1 _ 1 _ M 

^ ! Hwkrow(Sw) 1 + 1/M M+l 
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Diffusive flow - Buckley-Leverett theory 
Effect of gravity and flow rate 

► Gravity 

• Fractional flow equation shows that the gravity term will have a positive effect for oil 
displacement when 0 < a < n i.e. in the up dip direction and negative effect when 
n < a < 2n i.e. in the down dip direction 

• => the fractional flow of water for displacement in the up dip direction is lower than for 
displacement in the down dip direction since, in the first case, gravity tends to suppress 
the flow of water 

• => up dip is favorable to the displacement of oil by water 

► Rates 

• The fractional flow equation shows that the fractional flow f w decreases with 
decreasing rates 

=> low rates are favorable for displacement of oil by water 

=> slow production should be preferred 
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Diffusive flow - Buckley-Leverett theory 
Effect of capillary pressure 



The fractional flow equation shows that 
high capillary pressures lead to higher 
values of f w that are unfavorable to oil 
displacement 

However, at the front, high capillary 
forces tend to smooth the sharp 
solution resulting from Buckley-Leverett 
method 

=> The higher the flow rate, the closer 
the profile obtained by BL theory 
approaches the true profile 
=> In general, velocities in the reservoir 
are slow but the capillary forces can still 
be neglected 
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Diffusive flow - Buckley-Leverett theory 
Some ratio of interest 



► Fractional flow or water-cut 

q w 


fw = 


q w + q Q 1 + vJi row(^w) 

[L Q k 

rw C $w) 


= l-fo 


► Surface water-cut 


F w = 


:w 


q w /B 


w 


w 


Qw Qo qw/Bw “h qo/Bo ^ _|_ b_w ^ 

r / i 


► Water-Oil Ratio 


► Gas-Oil Ratio 


WOR = 


Qw _ q w /B w fw b 0 

Qo q 0 /B 0 i f w B w 

fg B g 


Qg 

Qo 


GOR = -^- = R s + 


1 ~fg B 0 
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Diffusive flow - Buckley-Leverett theory 



\ 

► Displacement of oil in the reservoir by water (or gas) is modeled through 
the assumption of so-called diffusive flow i.e. linear immiscible 
incompressible under vertical equilibrium 

► Buckley-Leverett theory provides a simplified solution for unidirectional 
displacement disregarding capillarity 

► Buckley-Leverett allows to calculate the fractional flow (also named water- 
cut) f w as a function of water saturation S w 

► Using the Welge tangent method, one can easily access to: 

• The estimate of the water saturation at the front S W f 

• The estimate of the mean water saturation behind the front S wm hence to the mean 
oil saturation after the front and corresponding oil recovery 

• The estimate of the velocity of the front of constant water saturation V f hence the 
time of water breakthrough 
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Sweep efficiency 
Definition 



► The sweep efficiency is defined as the recovery factor in reservoir conditions 
for the area that has been flooded: 


V s 

v p^oi 

where S oi is the oil saturation at the start of injection 


► The sweep efficiency can be expressed by: 

E = E a . E z . E d 

Where: 

E a is the areal sweep efficiency 

E z is the vertical (or invasion) efficiency 

E d is the displacement (or microscopic i.e. at the pore scale) efficiency 


► E is typically in the range 25-60% 
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Sweep efficiency 

Concept illustration for water flooding 
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Sweep efficiency 
Definitions 

► Displacement efficiency: 


E D = 



1—5 — s 

x ^wi ^ or 


► Vertical efficiency: 

vertical area swept 
total vertical area 

► Areal (superficial) efficiency: 

horizontal area swept 
a total horizontal area 





Production 
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Sweep efficiency 
Displacement efficiency 



► The displacement efficiency at the pore scale depends on: 

• Wettability 

• The capillarity effects and porous network 

• The mobility ratio: relative permeability & viscosity 



BP Videos - YouTube 

Exploiting science to increase oil recovery series 


In water-wet systems, at the pore scale, oil can be trapped in the center of the pores while water is 
flowing around. Oil that is connected to flow paths continues to be displaced. 

In oil-wet systems, oil may still be connected to flow paths but there may be large flow of water at the 
center of the pores. 
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Sweep efficiency 
Areal efficiency 



► The areal efficiency depends on: 

• The mobility ratio: relative permeability & viscosity 

• The pattern of injection 

• The permeability heterogeneities 

• The pressure distribution 



Oil can be bypassed because of the ineffectiveness of the macroscopic sweep. A pattern flood can be 
affected by a heterogeneous formation (especially anisotropic permeability i.e. preferential direction 
for permeability) or by the development of fingering/tongues because of a contrast of 
viscosity/mobility 

/PPTraining 
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Sweep efficiency 
Vertical efficiency 

► The vertical efficiency depends on: 

• The rock property variation between different flow units especially vertical distribution 
of permeability 



The vertical sweep can be affected by viscous fingering, as well as by the preferential movement of 
fluids along a high-permeability thief zone or by gravity override of injection gas or underride of 
injection water. 
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Sweep efficiency - Fingering 



► The Buckley-Leverett frontal advance theory assumes that the initial displacement of oil by water occurs 
as a smooth, substantially straight interface. 

► In 1951, Engelberts and Klinkenberg showed that, in scaled experiments, the existence of fingers, or 
discrete streamers of displacing water moving through the oil, could be inferred. Later van Meur 
demonstrated the existence of these fingers by a clever experimental technique. He showed that even in 
laboratory systems where care is taken to insure as nearly uniform porous media as possible, the 
tendency for these fingers to form increases as the oil-water viscosity ratio increases. 

► Later work by other investigators showed that at high oil-water viscosity ratios, there are instabilities at 
the oil-water interface and these instabilities grow until their effect is dominant on the overall flooding 
performance. 

► Debates raged over the practical significance of these fingers. An otherwise peaceful gathering of 
production research engineers could be turned into a heated controversy by the question, "Do you think 
viscous fingers are solely a laboratory phenomenon? 

► The actual reservoir having complex variations in permeability and porosity a macroscopically non 
uniform flood front very much similar to viscous fingering will occur, if from no other cause than the non 
uniformities in rock permeabilities. 

► A practicing reservoir engineer would term this the effect of reservoir heterogeneities. 
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Sweep efficiency 



5-spot areal efficiency 


fO 

OC 


-Q 

O 


O 

0 ) 

ro 

a> 

u 

Q) 

a 



producer 



► The superficial efficiency increases 
when : 

• the mobility ratio decreases 

• the velocity of the displacement 
decreases 


Increase of the Displacement Velocity 


Drive Mechanisms - Secondary Recovery 


From Petijeans 


!FP Training 


21 

6 


Sweep efficiency 

5-spot areal sweep efficiency and mobility ratio 



Sweep efficiency for the five-spot pattern as a 
function of volumes injected and mobility ratio 
M = Yw/Yo 

after Dyes, Caudle, and Erickson, Trans. AIME 



Reciprocal mobility ratio 1/M 


Sweep efficiency for the five-spot pattern as a 
function of fractional flow/ w and mobility ratio 
M = Yw/Yo 

after Dyes, Caudle, and Erickson, Trans. AIME 



Reciprocal mobility ratio 1/M 
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Sweep efficiency 

5-spot areal sweep efficiency and mobility ratio 




Mobility Ratio 


From Caudle 


At breakthrough (i.e. when the displacing fluid breaks at the well i.e. / = 0) it can be 
seen that E A = 0.7 for M = 1 while E A -> 0.4 for M > 1000 
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Sweep efficiency 

Effect of vertical permeability distribution 





Sweep efficiency 

Effects of vertical permeability distribution - 2 
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Sweep efficiency 



CXu 

\ 

► Sweep efficiency is defined as the recovery factor at reservoir conditions 
for the area that has been flooded 

► It is separated in three terms: 

• Microscopic or displacement efficiency related to the pore scale 

• Areal efficiency and vertical efficiency related to the flooding 

► Microscopic efficiency depends on the initial water saturation and the 
residual oil saturation 

► Areal efficiency depends on the mobility ratio and the water cut 

► Vertical efficiency depends mostly on vertical permeabilities and vertical 
distribution of permeabilities; due to heterogeneities, it is most difficult to 
calculate and is therefore very often estimated through reservoir 
simulation 

► Sweep efficiency is typically in the range 25-60% 

v 
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Generalized 
material balance 
equation for 
secondary 
recovery 


Generalized material balance 

Oil reservoirs generalized MBE with water/gas injection 



NB oi = initial oil volume 

(N — Np)B 0 oil volume left in the reservoir with its dissolved gas 

+ [NJ? S j — (N — Np)R s — G p ]B g volume of gas released from oil and staying in the reservoir 

+mN B oi {B g / B gi — l) volume of gas from the initial gas cap invading the oil zone 

water entry from the aquifer 
produced water 
injected water 
injected gas 
With R p = G p /N p the average produced GOR 

And B t = B a + ( R si — R s )B g the total or two-phase formation volume factor 


+ W, 


-W v B w 


+ W ■ B ■ ■ 


GinjB ginj 
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Generalized material balance 
Drive index 



Introducing R p within the previous equation and rearranging the terms: 


N V \B 0 + ( R p — R s )B g ] — N[(B 0 — B oi ) + ( R S i — R s )B g ] 

+mNB oi [B g /B g t - l] 

+w —WB + W B ■ ■ + G B ■ ■ 

~ vv e vv p^w ' vv injL'winj • ' J inj L, ginj 


► Introducing D = N p [B 0 + (R p - R s )B g ] = N p [B t + (R p - R si )B g ] and 
dividing by D we get: 

^ mN B oi[B g / B gi~ l] ^ W e ~WpB w ^inj^winj GinyZ?ginj 


D 

D 

D 

D 

-- DDI + 

SDI 

+ WDI + 

IDI 

Depletion 

Segregation 

Water 

Injection 

Drive Index 

Drive Index 

Drive Index 

Drive Index 
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Generalized material balance 
Injection drive index 

► By definition 

W B ■ ■ G B ■ ■ 

IDj _ vv iwy^winj gva] 

~ D D~ 



where 


D — N p [B 0 + (i? p — Rs)B g ] — N p [B t + — R S i)B g ] 


IDI allows to compare the relative importance of injecting water and/or gas as 
a drive mechanism by respect to the other natural drive mechanisms of the 
reservoir 




Secondary recovery 
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► Secondary recovery aims to improve oil recovery when natural drive 
mechanisms are not efficient enough i.e. when there is no strong natural 
water drive 

► It consists in supporting pressure and sweeping the remaining oil in the 
reservoir by injecting fluid, either water or gas 

► Secondary recovery implies multiphase flow in the reservoir that has to be 
controlled in order to get a high sweep efficiency 

► Sweep efficiency is the Rf at reservoir conditions and for the flooded 
area; it is divided into the microscopic efficiency, areal efficiency and 
vertical efficiency 

► Sweep efficiency depends on multiple parameters among which the 
mobility ratio M: if M > 1, it is unfavorable, if M < 1, it is favorable 

► Sweep efficiency is typically in the range 25-60% 
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Secondary Recovery 



<Xu 

I* * 

► 

► 


► 

► 


Feasibility studies are mandatory before performing a secondary recovery 

Economics will decide on the type of secondary recovery to be used 

• It is necessary to show an increase in the reserves associated with good economics 

• The impact on the surface facilities (pumping/compressing and treatment) has to be 
assessed 

The natural drive mechanisms have to be identified as soon as possible 

The distribution of heterogeneities in the reservoir, especially regarding 
permeability, have to be understood and modeled through reservoir 
modeling and simulation 

Extensive laboratory tests have to be done 

SCAL (wettability, Kr's, Pc) 

Water flooding and gas flooding at reservoir conditions 
Water compatibility issues 
Studying clay swelling issues 
Fine mobilization issues 
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Water injection 



\ 

► Water injection is the most commonly used secondary recovery process 
especially because of water availability 

► The mobility ratio may be favorable, especially for medium oil and it is, in 
any case, less unfavorable then for gas injection; however residual oil 
saturation may remain high 

• Density contrast is lower than for gas injection but high dip reservoir may favor 
gravity drainage and stabilization of the injection front 

► A good knowledge of the reservoir is mandatory: 

• Water/oil contact: injection should be done in the water zone for a better injectivity 
or at least near the contact 

• Permeability heterogeneities especially vertical permeability heterogeneities 

► The recovery factor can be as high as 50-60% 
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Gas injection 


0^ # \ 

► Gas injection is less used than water injection mainly because of gas 

availability and because of the cost of surface facilities (compressor) 

• The mobility ratio is always unfavorable for gas injection, instabilities develop quickly 
but oil residual saturation is lower after gas flooding than after water flooding, and a 
higher mobility of gas allows to flood areas remained unswept by water 

• Density contrast is higher and may favor gravity drainage and stabilization of the 
injection front 

• Immiscible gas injection can be turned to miscible gas injection with even lower 
residual oil saturation (possibly nil when miscibility is effectively reached) 

P A good knowledge of the reservoir is mandatory: 

• Gas/oil contact: injection should be done close to the contact if any 

• Permeability heterogeneities especially vertical permeability heterogeneities 

• A high dip reservoir may further favor gravity drainage 

P The recovery factor can be as high as 60-70% especially in case of gas 

gravity displacement 
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Recovery factors 


CXu 


Natural drainage mechanisms for oil reservoirs: 

• Monophasic expansion: Rf is a few % 

• Solution gas drive: Rf 10-25 % 

• Gas cap drive: Rf 25-40 % 

• Natural water drive: Rf 40-60% 

• Compaction drive: Rf 0-20% 


► Secondary recovery in oil reservoirs: 

• Water injection: Rf up to 40-60% 

• Gas injection (gravity displacement) : Rf up to 60-70% (see EOR) 
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Annex - Dry gas 
recycling in 
retrograde gas 
reservoirs 



Dry gas recycling 

Principle of dry gas recycling - 1 

► Retrograde condensate gas reservoir 

• In these reservoirs, when the pressure drops below dew point pressure, some liquid 
forms in the reservoir especially at the vicinity of the wells => phenomenon of 
condensate banking/blocking 

• Lower recovery for the liquid since a large part of it remains immobile in the reservoir 
due to a generally high value of critical condensate saturation 

• Lower recovery for the gas since the presence of the immobile liquid at the vicinity of 
the well disturbs the flow of gas 

► One solution: dry gas recycling 

• After separation at surface, all or part of the dry gas is reinjected in order to maintain 
pressure above the dew point 

• Reinjection of all of the dry gas: full recycling 

• Reinjection of part of the dry gas: partial recycling 

• Process is miscible but can still be modeled through black-oil modeling (oil & gas) if 
pressure does not fall below the dew point 
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Dry gas recycling 

Principle of dry gas recycling - 2 

► Full recycling 

• All the dry gas is reinjected but condensates are sold => deficiency in the material 
balance leading to a decline in pressure 

• Main economic parameter: initial condensate yield above the dew point r si 

• Recycling stops when the dry gas breaks at the well and the condensate production rate 
cannot sustain the required economics => blow-down and shift of sales towards dry gas 

• Dry gas injection may also limit water entries in case of an active aquifer thus limiting 
the quantity of gas trapped behind the front 

► Partial recycling 

• Typically reinjection of dry gas in excess of gas sales contract (e.g. during summer 
periods) 

• However a limited pressure maintenance may improve the recovery of condensates and 
gas and limit water entries if any 
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Dry gas recycling 
Mobility ratio 

► By definition 

_ krgdmax l^rgwmax 

► Relative permeability 

• Wet gas displacement by dry gas is a miscible process and IFT is nil 
=> S rgw = 0 in the pores contacted by the dry gas 

=> k r gd max — krgwmax 1 

► Viscosity 

• Pgd < Pgw but of the same order of magnitude typically 0.02 and 0.03 cP 

► As a consequence, M is always slightly higher than 1 but the process remains 
quite stable at macroscopic scale 
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Dry gas recycling 
Gravity effect 

► The table below shows the typical properties of dry and wet gas 



Wei gas 

Dry gas 

Gravity, (air = 1) 

0.937 

0.672 

Pseudo-critical temperature [2], l' c (PR) 

426 

374 

Pseudo-critical pressure [2], p c (psia) 

641 

663 

Z -factor at initial presure 

0.894 

0.936 

Gas expansion factor at p ly E\ (scfVrcf) 

267 

255 

Gas FVF at p\, B g \ (rcf/scf) 

0.00375 

0.00392 

Gas density in the reservoir (Ib/cu.ft) 

19.1 

13.1 


• FVF gd > FVF gw => this property allows to compensate for the effect of removing the 
condensate on the material balance 

• Pgd < Pgw => since M is slightly unfavorable, the stability of the displacement will 
strongly depend on the heterogeneities in the reservoir and especially the distribution of 
vertical permeability 

• A favorable case is coarsening downwards i.e. higher values of the permeability at the 
bottom 
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Dry gas recycling 
Vertical equilibrium 



► With no restriction of vertical fluid movement across the reservoir, 
segregation occurs with lighter dry gas rising up to the top of the section 


Coarsening upward of permeability will 
produce an early breakthrough of the 
dry gas and result in a poor sweeping 
efficiency => recovery should be low 
despite miscibility and large quantities of 
dry gas will have to injected 



Thickness 


fa) 


Coarsening downward of permeability 
will produce a very efficient piston-like 
displacement across the reservoir 


Permeability 



Cb) 


Injection 


Production 
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Dry gas recycling 
Pseudo-relative permeabilities 



► Despite miscibility, calculating pseudo-relative permeabilities in order to 
solve the displacement of the wet gas by the dry gas is still valid 

► In the case of N layers (ranging from top): 

• The average is done after the n first layers (from top) have been flooded 

• Each layer is assumed to have been completely flooded with a final dry gas saturation of 
1 — S wc . (i.e. residual saturation in wet gas is nil) or not to have been flooded at all 

• Thickness averaged water saturation is: 

s gdn = Z?=1 M ,(i - S W J/Y.U Mi 


• Thickness averaged relative permeability to dry gas is: 

^ rgd n ~ ^i=l^i^i^rgdmax/^i=l^i^i 


• Thickness averaged relative permeability to wet gas is: 

It — h.lr.lr /Y n h -It. 

n 'rgw 1 i ^ i=n+l ‘ L i n 'i n 'r gwmax / Zji=l 
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Dry gas recycling 
Vertical efficiency 



► Principle 






The generation of pseudo (averaged) relative permeabilities makes the displacement 
one-dimensional 

Buckley-Leverett theory and Welge tangent approach are suitable => the fractional flow 
is given by: 



1 ~ Hf fcr q W -0 sinO 


1 + 


A Igd. kygw 
Pgw k rgd 


1 - G 

1 _l_ K-gw 

Hgw k rgd 




Generally G can be neglected in the case of dry gas recycling 
Recovery: 

n _ $ gde “h (l — f gde)^id 

pD ~ w 

A J wc 

where G pD is the cumulative recovery of wet gas (in HCPV) 
and G iD is the cumulative dry gas injected (in PV) 
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Dry gas recycling 



0^ 

► Dry gas recycling aims to inject dry gas in retrograde gas condensates 
reservoir in order to maintain pressure 

► The objective is 

• to prevent the pressure to drop below dew point 

• to prevent the phenomenon of condensates banking in the vicinity of the wells and 
to improve recovery of condensates 

► Despite being a miscible process, the immiscible approach is still valid for 
dry gas recycling using black-oil model 

► Using averaged pseudo relative permeabilities and fractional flow makes 
the flow one-dimensional and allows to use Buckley-Leverett theory and 
Welge tangent method in order to forecast recovery 

► The final performance strongly depends on the vertical distribution of the 
permeabilities in the reservoir 

v ► The driving parameter is the liquid yield r 5 and the corresponding 
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